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Periodic solutions of planetary systems 
with satellites and the averaging method 
in systems with slow and fast variables 

The partial case of the planar N + 1 body problem, TV ^ 2, of the type 
of planetary system with satellites is studied. One of the bodies (the Sun) is 
assumed to be much heavier than the other bodies ("planets" and "satellites"), 
moreover the planets are much heavier than the satellites, and the "years" are 
much longer than the "months". Under a nondegeneracy condition, which in 
general holds, the existence of at least 2 N ~ 2 smooth 2-parameter families of 
symmetric periodic solutions in a rotating coordinate system is proved such 
that the distances between each planet and its satellites are much shorter than 
the distances between the Sun and the planets. Generating symmetric periodic 
solutions are described and necessity of the nondegeneracy condition is proved. 
Sufficient conditions for some periodic solutions to be orbitally stable in linear 
approximation are given. Via the averaging method, the results are extended 
to a class of Hamiltonian systems with slow and fast variables close to the 
systems of semidirect product type. 
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In memory of Nikolai Nikolaevich Nekhoroshev 
§ 1. Introduction 

We study the partial case of the planar N + 1 body problem, N 2, that can be 
characterized as "the problem on the motion of a planetary system with satellites". 

An effective estimate for the number of smooth two-parameter families of symmetric 
periodic solutions of this problem in a rotating coordinate system is proved (theorems 
1.1, 1.2(A) and corollary 1.1(3) about "solutions of the first kind"). Sufficient 
conditions for orbital stability in linear approximation for some of these solutions 
are given (theorem 1.2(B)). Generating symmetric periodic solutions are described 
(theorem 1.1). The necessity of a nondegeneracy condition is proved (theorem 1.3 
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and corollary 1.1 (:$)). The periodic solutions under our investigation are close to 
collections of independent "circular" solutions of the corresponding Kepler problems 
for each planet and each satellite. Via the averaging method (theorem 2.5 and 
corollary 2.1), the listed results are generalized to a wide class of Hamiltonian 
systems with slow and fast variables (theorems 2.1-2.4). 

Theorems 1.1 and 1.2 of the present work include as partial cases results of 
G. Hill [1, 2, 3] and H. Poincare [4] on the existence of periodic solutions and 
sufficient conditions of their orbital stability in linear approximation for the systems 
of the Sun-Earth-Moon and the Sun-two planets types (respectively). Theorem 1.1 
implies the known results by G. A. Krassinsky [5] and E.A. Kudryavtseva [6, 7] on 
the number of periodic solutions of planetary systems without and with satellites 
(respectively), by V.N. Tkhai [8] on the number and the location of symmetric 
periodic solutions of the systems of the Sun-planets and Sun-planet-satellites types. 

Let us formulate the results of the paper more precisely. 

A general N + 1 body problem is described by the system of equations 

where /i; is the mass of the ith body, is its radius vector in a fixed Euclidean 
space E, 

U = - (2) 



is the Newtonian potential of bodies' attraction, m = |r^ — r,-| (0 ^ i,j ^ N) are 
pairwise distances between the bodies, and g > is the gravitational constant. We 
consider a planar N + 1 body problem, i.e. the case of diraE = 2. 

Remark 1.1. Without loss of generality, we may choose the unities of mass, 
distance and time as it will be suitable. In fact, for any constants a,b,c > 0, a 
collection of vector-functions r,(t) is a solution of the N + 1 body problem with 
the gravitational constant g and the masses \ii if and only if the collection of 
vector-functions r,;(i) := b 2 cri(at/b 3 ) is a solution of the N + 1 body problem 
with the gravitational constant g = a 2 g and the masses jttj = c^. In particular, we 
do not need to assume that the gravitational constant g is arbitrary, but we may 
assume its value to be a distinguished number that we will choose below. (This can 
be achieved via scaling the time.) 

Definition 1.1. A solution of the planar N + 1 body problem will be called 
relatively periodic (or simply periodic) if the locations of all bodies after the time- interval 
T > can be obtained from their initial locations by rotating the plane by the same 
angle a around the centre of masses, for any initial time, where —it < a ^ n. The 
pair of real numbers (T, a) will be called the relative period of the solution, and the 
solution itself will be called (T , a) -periodic. 

Any solution obtained from a (T, a)-periodic solution via shifting the time by a 
value t and rotating by an angle ip around the origin is a (T, a)-periodic solution 
too. The union of the phase trajectories of all such solutions is a two-dimensional 
torus in the phase space, since it admits angular coordinates t modT, </?mod 2ir. All 
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these solutions will be regarded as a single (T, a) -periodic solution, and the union 
of their phase trajectories will be called the phase orbit of this solution. 

Many relatively periodic solutions of the planar N + 1 body problem happen to 
be symmetric (theorem 1.1). These solutions are characterized by the following 
property: at some time, all the bodies lie on the same line (i.e. a "parade" is 
observed) and their velocities are perpendicular to this line. 

In the present work, the following partial case of the planar N + 1 body problem 
is considered, N ^ 2. We assume that the mass of one particle (the Sun) equals 
/io = 1 and is much greater than the masses of the other particles (the planets 
and satellites) . Moreover the mass fa of the ith planet and the mass faj of its jth 
satellite have the form 

fa = firrii, faj = [ivm, l3 < fa, U i ^ n, 1 < j < fa, (3) 

where < fa v -C 1 are small parameters and mi,mij are positive parameters far 
enough from zero (e.g. positive constants) with the properties 

n Tii 

2> = 1, min£^ = L (4) 

z=i j=i 

moreover my are bounded for fa ^ 2, where rii is the number of satellites of the 
ith planet and 1 + n + Y27=i n>i = N +1 is the number of all bodies. Thus, for each 
"double planet" (rii — 1), the mass of the satellite equals fiiriidi/(l — Qi) where the 
parameter 6i :— vmn/(mi + vmn) € (0, 1/2] is not necessarily small (since ran is 
not necessarily bounded). Due to (4), one has 9i(l — 6i)jv = (1 — B i)mn j fhi ^1/4 
for rii = 1. 

We also assume that the distance Ri between the Sun and the ith planet is 
of order R 3> 1, while each satellite is at the distance ry of order 1 from its 
planet. Finally, "the years are much longer than the months", i.e. the angular 
frequency uii of the rotation of each planet around the Sun is of order w, while the 
angular frequencies fiy of the rotations of its satellites about it have order 1 where 
< oj -C 1. More precisely, let a set of non- vanishing real numbers 

uji = u)fl i0 , flij, K K n, l^j^m, (5) 

called the set of angular frequencies satisfies the conditions 

c< |fi i0 | < \SUj\ < 1, l^i^n, l^j^fa, (6) 

||fi i0 | - l fi i'o|l > c, i< i', \\tUj\ - 1%'H > c, 1 s$ j < f sC m. (7) 
Here c is a suitable real number in the interval < c < 1. 

Suppose that the parameters u,^,,g,R satisfy the natural relations w 2 i? 3 = g 
and 1 = gfj, corresponding to Kepler's second law for the planets (wfi?f = g^o) 
and the satellites (fi^-rf^ — gfa), for the chosen unities of mass, distance and time. 
Thus, g = w 2 R 3 = 1/fa 

^3=^V (8) 

and the problem has three independent small parameters: fa v and w. We emphasize 
that the initial N + 1 body problem does not include the parameters to and R, so we 
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introduce them as additional small and big parameters. Therefore the only imposed 
restrictions to the parameters of the problem are as follows: the mass of the Sun 
is fio = 1, the distances from the satellites to their planets are of order 1, and the 
gravitational constant is g = 1///. This does not cause any loss of generality due to 
remark 1.1. 

Let us describe the "generating" periodic solutions of the N + 1 body problem 
under consideration. These are the periodic solutions of the "model" problem (see 
remark 3.2 and the formula (73)) consisting of N independent Kepler's problems, 
one for each planet or satellite. Let us assume that the collection (5) is maximally 
resonant, i.e. has the form 

% = uii+Kij^-, 1 <i<ni, 1 ' 

where ki,Kij £ Z, T > 0. The solutions of the model problem will be called 
generating solutions if they correspond to independent circular motions of the 
planets around the Sun (which is placed to the origin) and the satellites around 
the planets with the angular frequencies (5) of the form (6), (7), (9). Here the 

1/3 

circular orbits of the planets and satellites have radii Ri = ^ ^j 2/3 and = ^"^2/3 
respectively, 1 ^ i ^ n, 1 ^ j ^ n^. The union of the phase trajectories of all 
such solutions is an iV-dimensional torus, since the polar angles of N radius vectors 
drawn from the Sun to the planets and from the planets to their satellites can be 
used as coordinates on it. Denote this torus by A°. Due to the condition (9), 
generating solutions are (T, a)-periodic with relative period 

T>0, a = Lo 1 T + 2irk £ (-7T,7r] (10) 

where k £ Z is a suitable integer. In the presence of satellites, the following condition 
will be called the nondegeneracy condition: 

\a\>io 2 T. (11) 

Symmetric periodic solutions of the three-body problem were studied already by 
Poincare [4] . Recall the definition of a symmetric solution of the planar N + 1 body 
problem. 

Definition 1.2. Consider a problem describing the motion of N + 1 particles 
in a Euclidean plane. A solution of this problem will be called symmetric if there 
exists a line I in the plane, called the axis of symmetry, and a time t — to satisfying 
one of the following (equivalent) conditions called a "parade" of the particles: 

1) at the time t — to, all points are placed on the line I (i.e. a "parade" of the 
particles is observed) and their velocities are orthogonal to the line 

2) the locations (and, hence, also the velocities) of all particles at any time t £ R 
can be obtained from their locations at the time 2to — t by reflecting with respect 
to the axis I. 

The particles of the system are assumed to be numbered. The order of them on 
the line I at the time to can be called the type of a parade. Any solution of the 
N + 1 body problem obtained from a symmetric solution by shifting the time and 
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by rotating the plane is also symmetric. Similarly to the case of (T, a)-periodic 
solutions, we will not distinguish such solutions and will regard them as a single 
symmetric solution. 

One easily shows that exactly 2 N ~ 2 of the generating solutions are symmetric. 
Namely, the symmetric (T, a)-periodic solutions are characterized by the condition 
that "parades" of the planets and satellites are observed, moreover they repeat each 
half of the period, T/2. This means that all the particles of the system are posed 
on a line, which turns by the angle a/2 after the time-interval T/2. 

Theorem 1.1. There exist constants ujq,C > and continuous positive functions 
fiQ = /io(w,c) and = i/o(cj,c) (0 < w ^ ujq, < c < 1) such that, for 
any parameter values u>,c,/j,,v with the properties < u> ^ ujq, < c < 1, 
< /! ^ /io( w i c ) an d < v ^ vq[uj,c) and any collection of angular frequencies 
(5) of the form (6), (7) and (9), the following property holds. Suppose that the 
parameters (10) satisfy either the nondegeneracy condition (11) or the following 
more delicate conditions: 

UJ 2 

a^O, a- — j-T # [-Cuj 3 T, Cuj 3 T] + 2ttZ (12) 

4:1 L<lj 

for 1 < £ < ri, 1 < j < rii- Then the N + 1 body problem of the type of 
planetary system with satellites, N ^ 2, has exactly 2 N ~ 2 symmetric (T, a) -periodic 
solutions that are O(oj 2 )-close to the generating symmetric solutions corresponding 
to independent circular rotations of the planets and satellites with the angular 
frequencies (5). Each of these 2 N ~ 2 solutions smoothly depends on the pair (T, a). 
For each of these solutions, parades are observed, which repeat each time-interval 
-j: all of the particles of the system are posed on a line (which turns by the angle 
a/2 after the time-interval T/2). 

Let A° C T*Q be the TV-dimensional torus in the phase space formed by the 
trajectories of the generating solutions (see (10)). Let S C T*Q be a "transversal 
surface" of codimension 2 in the phase space, which transversally intersects invariant 
two-dimensional tori lying on A and corresponding to (T, a)-periodic solutions: 

^2(<Pi + P«) = Yl^Wi +^2Kijcpij) = 0(mod27r) 

i=\ j=l i=\ j=l 

Here ki, Kij are integers in (9), T m j n is the minimal positive period, hence the integer 
T/T m i n is the greater common divisor of the collection of integers ki, K^. Let A be 
the A^-dimensional torus that is w 2 -close to the torus A° (see theorem 1.2 below). 
Let ^> be the generating function of the "succession map" of the N + 1 body 

problem under consideration (see definition 2.1). Consider the smooth function 
S = ^Ixns on t ne (N — 2)-dimensional torus AflE. Since the function S is defined 
on a (N — 2)-dimensional torus, it has at least N — 1 critical points, moreover at 
least 2 N ~ 2 points counted with multiplicities [9]. We will prove (see theorem 1.2) 
the same lower bound for the number of (T, a)-periodic solutions of the problem 
under consideration. Observe that the function S has at least one critical point, 
since it is defined on a closed manifold. We will prove that each critical point of 
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the function S corresponds to a (T, a)-periodic solution of the N + 1 body problem. 
Moreover, we will offer sufficient conditions that guarantee the orbital (structural) 
stability in linear approximation (see definition 3.1) of such a solution. 

The following condition will be called the strong nondegeneracy condition: 

u 2 T <\a\<n-uj 2 T. (13) 

The following conditions will be called the property of having fixed sign: 

1) all planets rotate "to the same side", i.e. the angular frequencies of the rotations 
of the planets around the Sun have the same sign: 

WjW.,' > 0, 1 < i < i' ^ n; (14) 

2) the function S on the (N — 2)-dimensional torus is either a Morse function or 
has at least one nondegenerate critical point of a local minimum (this condition is 
assumed to be always true if N = 2). 

THEOREM 1.2. Under the hypothesis of theorem 1.1, there exists a smooth N -dimen- 
sional torus A in the phase space that is uJ 2 -close to the torus A , smoothly depends 
on the pair (T, a) and has the following properties. 

(A) The phase orbits of all (T, a) -periodic solutions oftheN+1 body problem that 
are u-close to the torus A° are contained in the torus A. Moreover their intersection 
points with the transversal surface £ coincide with critical points of the function 
S := ^l^ns on the (N — 2) -dimensional torus AflE. Here $ is the generating 
function of the "succession map" gJ I _ UJl i of the problem under consideration (see 
definition 2.1). The function is an even function in the collection p of angle 
variables Pi\r, Pij\]{- The phase orbits of the symmetric (T, a) -periodic solutions 
contain the points if of the torus A having the property p = — ip. 

(B) Suppose that the property of having fixed sign holds, moreover either the 
strong nondegeneracy condition (13) or the following more delicate conditions hold: 



sgn ujx + sgn Q,ij 



-<L U »T ,-u 3 T 
2 ' 2 



sgnfiy +sgnfiiv / uf uf, U . , , 



ttZ, (15) 



a -{W\ + ]^)s^~ c ^ c ^ + ^ (16) 

for all 1 ^ ^ n, 1 ^ j ^ Ui and 1 ^ j' ^ ny . Then the (T , a) -periodic solution 
corresponding to any nondegenerate critical point of a local minimum of the function 
S is orbitally structurally stable in linear approximation. 

Theorem 1.2(B) implies that, for N = 3, in the "generic case", a half of the 
(T, a)-periodic solutions that are close to the torus A° are orbitally stable in linear 
approximation (since the function S is defined on a circle and, hence, has only 
critical points of local minima and maxima, which alternate on the circle) . 

The inequality |a| ^ it and the nondegeneracy condition (11) (respectively 
(12)) imply that, in the presence of satellites, the period T of the solutions under 
consideration is "not too big": respectively 
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Remark 1.2. (A) One easily checks that the nondegeneracy condition (11) is 
always realizable (for suitable collections of the angular frequencies u)i, £lij). In fact, 
for any < ui u>o(a) ■= y-, the period T can take any value of the form 

T = — , a > a (k 2 , ...,k n ):= max {7|fc;|, ^7(N - n + 1)} , (18) 

thus the value ui 2 T can be arbitrary small, hence any of the nondegeneracy conditions 

(11) and (12) can be always achieved. In more detail, let us fix any integers 
k%, • • • ik n € Z \ {0} having different absolute values, as well as any values a 
a (k 2 , . . . ,k n ) and c <E (0,c o (a)] where c (a) := min{^, 14(JV ^„ +1 ) }■ For any 
be [f , f] and cj e (0,w o (a)] put 

oji := u;b, uii := ui\ + , := cji + Kij— (19) 

for any integers -K^ G Z with the properties |l£y| € [|^, |^] and — 5* ^ 

(for all i, j 7^ /) where € g [^j, 7 ^ JV _° +1 ^ ] fl N. Such integers I and l£y exist 
for < uj < w (o) = since the interval [77, 7(]V ,° +1) J (respectively [%,%}) 

is of the length g( 7(Af _ 1 n+1) - c) > 7( jv-n+i) > 2 (respectively ^ ^ (AT - n + 
1)£) and, hence, it always contains a positive integer (respectively N — n different 
positive integers with pairwise distances £). Since a = uj\T = 2irab( mod 27r) 
and lo 2 T = 2imui ^ ^, the required nondegeneracy condition (11) has the form 
(0(6 - uj),a(b + uj)) n Z = 0. Hence it holds if (a& - |,a& + \) n Z = 0, i.e. 
6 € i + [-i, i] + iZ. It follows from ± > i that, for any a ^ a (fc 2 , . • . , k n ), 
there exists 6 € [1, |] PI ([j^, y-] + -Z). Clearly the nondegeneracy condition (11) 
holds for any such b. The obtained real numbers Cc>j,Qjj satisfy the inequalities 
c < i < M < I < |%| < 1, I -I"" 11 ^ I ^ c, HOyl - > |0 c and, 

hence, the inequalities (6) and (7). This proves the realizability of any period T of 
the form (18). In particular, for any 

a ^ a Q (k 2 , . . . , k n ), b = b(a):=- 

a 

the condition b £ [|, |] and the nondegeneracy condition (11) automatically hold for 
the period (18) and the collection of angular frequencies (19). Thus the system of 
relations (18), (19) and (20) (called the rough nondegeneracy condition) implies the 
nondegeneracy. Due to (6), the nondegeneracy (11) implies the delicate nondegeneracy 

(12) for small enough 0<uCl. 

(B) In the case n > 1 (of systems with more than one planet), the relative period 
T is always "big", namely T = , 27r|fc21 > 2JM ^ m > 1 for < u < 1, due to 
(9), (6) and (7). Together with any of the inequalities (17) (which are corollaries 
of the nondegeneracy condition (11) or (12)), this implies that the period T under 
consideration belongs to the interval (■^,■^7) or (^, -^jj) respectively. In the case 
N = 2, n = 1 (of a system of the Sun-Earth-Moon type), the minimal positive 
relative period is bounded and equals T m j n = 2n/\fln — u>i\ ^ Ait/c, while the 
rotation angle a m i n = wiT m ; n is of order uj. Hence the nondegeneracy condition 
(11) automatically holds for T = T m i u . In the case n = 1, N > 2 (of a system of the 
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Sun-planet-satellites type), the minimal positive relative period can be bounded 
too. However we do not assume in theorems 1.1—1.3 that the period T > in (10) 
is minimal. We only assume that it satisfies either the nondegeneracy condition (11) 
or the delicate condition (12) or the rough condition (18), (19), (20). In particular, 
we can assume that T = £2SSi ) see (A). 

The natural question arises: is the nondegeneracy condition (11) necessary for 
the validity of theorem 1.1? An answer happens to be affirmative in many cases. 

In the following theorem, by "almost any" collection of masses fa > 0, 1 ^ i ^ n, 
we mean any collection belonging to the complement in R™ to the union of a finite 
set of linear subspaces of R™. Moreover each of these subspaces depends on the 
collection of integers fcj in (9), has codimension at least 2, and the number of these 
subspaces does not exceed 2" -2 . The set A4 sym of "almost all" collections of masses 
is described in more detail in §4.1. By the phase space of satellites, we regard the 
direct product of big balls in the phase spaces of the corresponding Kepler problems, 
except for a small neighbourhood of "the set of possible collisions". 

Theorem 1.3. Suppose that, under the hypothesis of theorem 1.1, the number n 
of planets is at least 2 and there exist two planets whose angular frequencies satisfy 
the following resonance relation: 

fceZ\{-i,o}|. (21) 

In this case, a — automatically. Then there exist an open dense subset M Byva C 
R" and a nonempty open subset M C Ji4 sym (see definition 4.1 and remark 
4.1), both invariant under multiplication by any positive real number and having 
the following properties. For "almost any" collection of planets' masses fimi > 0, 
namely for any collection of planets' masses fi(mi, . . . ,m n ) G A^ sym (respectively 
for any collection of planets' masses (J,(mi, . . . , m n ) G M., for example satisfying 
the inequality |/«M'|c Kii ,mj' > v \ K u\ c nu' m h see (81), (82), (83)), there exist 

numbers (ig, vq > and an open subset Uq in the phase space of the planets containing 
the phase orbits of all symmetric "circular" solutions (respectively all circular solutions) 
of the collection of the Kepler problems for planets with angular frequencies (5), such 
that the following condition holds. For any values fx, u, T, a G R of the form 

°<(^) ^<Mo, \T-T\ + \a\^Dn, 

there exists no (T, a)-periodic solution of the N+l body problem under consideration 
whose phase orbit has a nonempty intersection with the direct product U of Uq and 
the phase space of the satellites. 

In particular, the region U does not contain the phase orbit of any symmetric 
T -periodic solution (respectively T-periodic solution). 

Consider the planetary system with two planets, a partial case of the three-body 
problem. In this case, the minimal positive period T m ; n equals i tl)2 _! r fa , 1 | > thus the 
condition (21) means that the corresponding angle a m ; n = — 27rfc vanishes. 

We also observe that, in this case, the region U in theorem 1.3 contains the whole 
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Fig. 1. A set of pairs (k = ^ , u) such that there exists (3) or 
does not exist ($) a relatively-periodic solution of the 3-body problem 

two-dimensional torus A° and M syra = M = M> (i.e. "almost any" means "any"). 
Thus, theorems 1.1 and 1.3 for N = n = 2 imply the following. 

COROLLARY 1.1. Consider the three-body problem of the type of planetary system 
with two planets. Fix angular frequencies of planets CJi ^ 0, W2 ^ 0, |a>i| =/= |k>2|> 
and consider the two-dimensional torus A° corresponding to the circular motions of 
planets with frequencies u>i,u>2- Put T := y^^^r, a. — 2ir + ^Trk £ (— 7r,7r] 
for a suitable k 6 Z. In dependence on the ratio of these frequencies, one of the 
following statements holds. 

(3) Suppose that the angular frequencies U)\,U)2 does not satisfy the special resonance 
condition (21). Then a/0 and there exists a number fiQ = no(rni,m2,oJi,u>2) > 
such that, for any values fi,uji,0J2, ImI + l^i — w i| + 1^2 — w 2| ^ Mo? ^ere exists 
a two-dimensional torus A^ 1 ^ 1 ^ 2 that smoothly depends on the triple (/i, Ci,^), 
coincides with the torus A° if (yU,w 1,^2) = (0,Wi,u;2) and has the following property. 
IfO<fi^fio then the torus A /ii 2i,S 2 * s ^ e phase orbit of a symmetric (T, a) -periodic 
solution of the problem under consideration with parameters T — ^-fi^ > ® = 
27r - — h 27rfc, for a suitable k S Z. 

(J) Suppose that the angular frequencies u>i,UJ2 are in a special resonance (21). 
Then a = and, /or any numbers T, D > with ^ 2 2 ^" 1 ^ G 2, t/iere eizsi a number 
Ho = jj,o(mi,m2,oJi,u}2,T, D) > and a neighbourhood U — U mijm2tU;itU2t T t D of 
the torus A in the phase space such that, for any parameter value /i G (0, /Iq] of 
the three-body problem under consideration, U does not contain any (T, 5) -periodic 
orbit with parameters T , a of the form 

\T-T\ + \a\^D(i. 

The figure shows (for a fixed u>i ^ 0) regions in the plane R 2 consisting of pairs 
(k, /j.) £ lx M>o, K '■= ul !f_l UJl 1 < /i ^ /io(mi,m2,wi,W2), such that there exists 
(respectively does not exist) a relatively-periodic solution of the three-body problem 
with parameters T = t^~^-| > a = 27tk + 27rfc, for a suitable k € Z. 

Due to corollary 1.1, for the planetary system with two planets, the condition 
(21) is false if and only if there exists a (T, a)-periodic solution close to a "circular" 
(T, a)-periodic motion. These (T, a)-periodic solutions were discovered already by 
H. Poincare [4] who called them solutions of the first kind. In the degenerate case 
(21), Poincare discovered periodic solutions corresponding to elliptic motions of the 
planets, solutions of the second kind. 

Theorems 1.1—1.3 seem to admit an extension to the cases of the N + 1 body 
problem on a sphere or on the Lobachevskiy plane, due to periodicity of solutions 
of the Kepler problems on these surfaces (see [10] and references therein). 
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The paper is organized as follows. In §2, theorems 2.1-2.4 on periodic solutions 
of Hamiltonian systems with slow and fast variables are formulated, and a scheme 
of their proof is sketched via generalizing the averaging method (see theorem 2.5 
and corollary 2.1) and the method of generating function (see definition 2.1). In 
§3.1 and §3.2, the parameters and the phase space of the N + 1 body problem under 
consideration are described in more detail and the notions of a symmetric solution 
and an orbitally stable in linear approximation solution are discussed (see definition 
3.1). In §3.3 and §3.4, we introduce coordinates in the phase space of the N + 1 
body problem, which bring the system to the form of the systems studied in §2 
(lemmata 3.1 and 3.3), moreover we establish relations between our unperturbed 
system, the Hill problem, the three-body problem and the restricted three-body 
problem (§3.3.2). In §4, theorems 1.1—1.3 are derived from theorems 2.1-2.4. 

A proof of the technical theorem 2.5, its corollary 2.1 and deriving from them 
theorems 2.1-2.4 will be published in our future work. 

As a conclusion, we remark that the question of interpreting the discovered class 
of relatively-periodic solutions of the N+l body problem in terms of the behaviour of 
planets and satellites of the real solar system is very interesting, needs an additional 
investigation, and it is not discussed in the present paper. 

The author is grateful to N. N. Nekhoroshev for stating the problem, A. D. 
Bruno, A. T. Fomenko, V. V. Kozlov, A. I. Neishtadt, V. N. Tkhai, Yu. M. Vorobiev, 
and H. Zieschang for their attention to this work, A. B. Kudryavtsev for numeric 
computation of the numbers c K from definition 4.1 of the sets M. C Ai sym C K™ 
in theorem 1.3. 

§ 2. The averaging method for a class 
of systems with slow and fast variables 

In §3 and §4, we will derive theorems 1.1-1.3 from the next theorems 2.1-2.4 on 
periodic solutions of dynamical systems having the following special form. 

Let p : M — > Mq be a locally trivial fibre bundle of smooth manifolds. Fix a 
number Ael. A pair of differential fc-forms (uj, Wi) on M will be called a X-pair if 
the fc-form u>q := u: — Xu>i "projects" to the base Mq of this fibre bundle (i.e. has the 
form ujq — p*u>o for some /c-form ujq on Mo), k ^ 0. A A-pair (o>, Wi) of 2-forms on 
M will be called a X-symplectic structure if the 2-form u)q is a symplectic structure 
on Mq, oji is closed and 

T X M = (Kcr 

at any point x € M. (Thus the 2-form u)\ determines a symplectic structure on 
each fibre, moreover its field of kernels is transversal to the fibre at any point of M 
and determines a "symplectic" flat connection of the fibre bundle p) . In this case, 
the fibre bundle p will be called a symplectic semidirect product and denoted by 
(M, Mq,p] u>0; k>i)- A vector field v = vh.Hi 011 M will be called X- Hamiltonian if 

(U (;V) - dH) iKor^i = 0, («].(•, - dH{) | K erw = 



for some A-symplectic structure (u>,u>i) and A-pair of functions (H,Hi) on M. In 
this case, the dynamical system x(t) = v(x(t)) on M will be called X-Hamiltonian 
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and denoted by 

(M,Mb,p; « ,a;i;ff,ffi) A , (22) 

moreover the functions H and Hi will be called the Hamiltonian function and the 
X-Hamiltonian function (respectively) of the system (22). If A 7^ then the system 
(22) is equivalent to the Hamiltonian system (M,uj,H): 

(M,Af ,p;u>o,«i;lf,ifi) A = (M,e>,H). 

Here the symbol = denotes the equivalence of (A-)Hamiltonian systems, i.e. the 
coincidence of the corresponding (A-)Hamiltonian vector fields. If A = then 
the system (22) is a semidirect product. (Such systems are studied by Yu. M. 
Vorobiev [11]; they include e.g. the restricted three-body problem, see §3.3.2). 

Denote by g^j H the flow of the A-Hamiltonian vector field Vh,Hi ■ Similarly to 
the case of Hamiltonian systems, the flow of the field Vh,Hi always preserves the 
2-form us and the Hamiltonian function H. 

As a very special example, consider the case when M = Mq x M% is a direct 
product, moreover the 2-form u}\ and the function H\ "project" to Ml (i.e. have the 
form u>i — p\uj\ and Hi = Hi o pi for some 2-form <2>i on M\ and function H± G 
C°°(Mi) where pi : M — > Mi is the projection). Then the A-Hamiltonian system 
(22) for any A € R is equivalent to the Hamiltonian system (M, u>q + lji,Hq + Hi), 
i.e. to the direct productofthe Hamiltonian systems [Mq, u>q, Hq) and [M\, cDi, Hi): 

(M, M ,p; u)q, U\\ H, Hi)^ = (A/, wo + u>i, Hq + Hi). (23) 

Let us describe a class of dynamical systems with slow and fast variables. Consider 
a two-parameter family of e-Hamiltonian systems on M with £-symplectic structure 
(u>, u>i), the Hamiltonian H = luHq+sHi and the £-Hamiltonian Hi where |e|, \w\ <C 
1. Here the 2-forms u»o,Wi and the functions Hq = Hq op, Hi depend, in general, 
on the small parameters e, ui and possibly on some other parameters of the system, 
moreover some relations between parameters may be posed. The local coordinates 
of a point Xq G Mq are "slow variables", while the local coordinates of a point "on a 
fibre" are "fast variables" of the system. 

Example 2.1. The following problems are described by systems with slow and 
fast variables: 

(i) the problem on the motion of a charged particle in a strong magnetic field on 
a symplectic manifold (Mo,u>o) where M — T*Mq is the cotangent bundle, the 
magnetic field if given by the 2-form e _1 / 4 o)o on Mq, and the small parameters lu, e 
are related by the condition ui 2 = e; 

(ii) the N + 1 body problem of the type of "planetary system with satellites" where 
M = Mq x Mi is the direct product of the phase spaces of planets and satellites, 
uji = p*Qi, lo = u>q + slji, and the Hamiltonian function is 

H = ljHq(xq;h) + e (Hi{xi\v) + uj 2 <S>(xq, xi, n, v, p)j . 

Here (xq, xi) S Mq x Mi, < ui, e, \i, v, p <C 1 are small parameters related by the 
conditions p = u; 2 / 3 /// 3 and e = u; 1 / 3 /^ 2 / 3 ^. The planar N + 1 body problem is 
^-symmetric and reversible (see below). 
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From now on, we assume that M = Mo x Mi is the direct product and u>i = p\(^i ■ 

Suppose that each symplectic manifold (Mi, <Dj) is equipped with the Hamiltonian 
action of a circle 50(2) = 5 1 = R/27rZ with the Hamiltonian function 7j, i = 0, 1. 
The system (22) will be called S 1 -symmetric (or SO (2) -symmetric) if the functions 
H, Hi are invariant under the diagonal action of the circle on M (i.e. invariant under 
the flow of the A-Hamiltonian field vjj 1 on M where I = Iq + Aii). All solutions of 
this system that differ by shifts along the commuting vector fields vh,H! and vjj x 
will be regarded as a single solution, and the union of their phase trajectories will 
be called the phase orbit of this solution. Let T, a £ R, T ^ 0. A solution 7(4) of 
a $ -symmetric system will be called (T, a) -periodic if it is defined on the whole 
time-axis, and j(t) = ffg jj^J" (l(t)) for some (and, hence, for any) teK. 

A ^-symmetric system (22) will be called reversible (or O (2) -symmetric) if each 
Mi is equipped with an anti-canonical involution Ji : Mi — > Mi preserving the 
function Ii (i.e. J*u>i = — u>i and IioJi = Zj), i = 0, 1, moreover the functions H, H\ 
are invariant under the (component-wise anti-canonical) involution J := Jo x J\ : 
M — » M. A solution j(t) of the reversible system will be called symmetric if it is 
defined on a time- interval (to — e, to + e) C K and 7 (to) = 9j°i 1 J gj^° (7(^0)) for some 
ipo € Rmod27r (and, hence, 7(2i — t) = .9/'} Jffj j°(7(t)) for any t 6 (t — io+s))- 

Let us describe the model system on M = Mo x Mi: 

(M, M , p; w 0) wi ; uH , H x f « (M, w + w 1} wff + H x ), (24) 

< cj < 1 where Jfi = #1 opi for some function H x € C°°(Mi), see (22) and (23). 
We will assume that each of the Hamiltonian systems (Mo, cDo, -ffo) and (Mi, u>i, Hi) 
in the system (24) is the direct product of S^-symmetric Hamiltonian systems: 

n n ni 

(M ,Qo,H )=Y[(M i0 ,Q i0 ,H i0 ), (M x , Q x , H x ) = JJ \{(M l3 , u>y, H l} ). 

i—l 

Moreover each factor My = S 1 x (ay,6y) x R 2 is equipped with coordinates 
ifij mod2n,Iij,qij,pij such that tDy = <iJy A cft^y + dpy A day, the action of the 
circle on (My,Wy) is given by the Hamiltonian Zy, and the involution J acts 
component-wise in the form (</?y , Iy , qy ,f>y) H- (— yy , Zy , ay, — Pij)- In particular, 

n n ni 

# = 51 #i0 ) -ff 1 = ^ ^ , TZy = Hij (Zy , Oy , py ) 

2=1 Z=l j=l 

where, from now on, we use the same notation for a function and its lift by misuse of 
language. We will assume that each coordinate cylinder 5* 1 x (ay, fey) x {(0,0)} C 
My consists of relative equilibrium points, i.e. the co- vector iffy is proportional to 
dlij at any its point (with coefficient depending on the point), 1 ^ i n, ^ j ^ rii. 
Therefore the 2iV-dimensional symplectic submanifold 

n rii 

J] jj S 1 x (ay , bij ) x {(0, 0)} C M x Mi (25) 

i=ij=a 

is invariant under the flow of the model system (24) and it is fibred by invariant 
iV-dimensional tori n"L S 1 x {(Zy, 0, 0)} where N := £" =i (l + ni). Those 
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solutions of the model system whose phase orbits are contained in the invariant 
submanifold (25) will be called the generating solutions. Consider one of these 
TV-dimensional tori, A°, and a (4./V — 2)-dimensional "transversal surface" S in the 
4iV-dimensional phase space Mq x Mi, which is transversal to the two-dimensional 
phase orbits of the generating solutions contained in the iV-torus A° . 

Let us describe the unperturbed system. Suppose that a ^-invariant function 
Fij = Fij(I i0 ,qi ,p i0 , ipij—ipi , kjiQijiPij) is given on each direct product M io xMy, 

l ^ j ^ m. Put 

* : =EE^- ( 26 ) 

1=1 j=l 

As the unperturbed system, we will regard the O-Hamiltonian system 

(M x M^Mo.p^o^^wffo.ffi + w 2 *) (27) 

with parameter < co <C 1- Then the S^-action is given via the Hamiltonian 
function Iq = I<o an d the O-Hamiltonian function l\ = X)"=i 

Both systems described above: the model one (24) and the unperturbed one 
(27), are systems with slow and fast variables, because of a small factor u> in their 
Hamiltonian function luHq. 

Theorem 2.1 (the number of relatively-periodic solutions). Suppose 
that all solutions of each Hamiltonian system (My, cDy, Hy) are periodic with periods 
Tij o H^, for some functions 7y — Tij(h) ^ 0. Moreover the number of the systems 
is N ^ 2. Suppose that every function H^ — Hij(Iij,qij,pij) satisfies the following 
conditions at all points (7^, 0,0) G (ciij,bij) x {(0,0)}: 

dHi^iUjVijWy, ^|^0> det g£^y 2 =«?j(Jtf) (28) 

where flij(Iij) := 2n / T t j (Hi 0,0)), 1 < i < n, ^ j ^ n t . 

Then, for any collection of segments [a°j, b°j] C (oy, bij), there exist real numbers 
luq,Ci,C2 > such that the following conditions hold for any ui € (0,Wo]- Suppose 
that, for some numbers 1°^ S \a°^,b°^, the collection of 11 angular frequencies" fiy = 
fly(J?-) satisfies the following two conditions. The first condition is the "relative 
resonance" condition: 

— - 2ir 

u) ~ J Qij = u>i + hj — , 1 < i < n, O^j^ni, (29) 

where integers kij £ 1 do not vanish simultaneously, t :— max{0,£}, uj\ := wfiio 
and T > 0. The second condition is either the nondegeneracy condition: 

a := wiT [~duj 2 T, du 2 T] + 2ttZ, (30) 

or the following more delicate condition: 

a 2ttZ, a + A^uj 2 T [-C* 2 u; 3 T, C 2 oj 3 T] + 2nZ (31) 
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for all 1 ^ i ^ n, 1 $C j ^ m. Here the real number Ajj is expressed via the first and 
second partial derivatives of the function (F°j) = (F°j} {qij , Pij) at the point (0,0), 
for instance 

Here the function (F°j) = (F°j) (Qij > Pij) * s obtained by averaging the function F°j := 

o 0)}xJ?r. 1 (H- •(/?. o o)) a ^ on 9 ff~ -periodic solutions of the system (Mij, Qij, Hij), 
C\ > |Ajj|, 1 ^ i ^ ii, 1 ^ j ^ Uj. T/ien £/iere exists £o > smc/i £/iai, /or 

< e ^ £o, any S 1 -symmetric ("perturbed") Hamiltonian system 

(M, w + eu!,H) = (M, M ,p; w„, w x ; 5, #i + w 2 $) £ (33) 

on M = Mo x Mi /ias at least N — 1 (T, a) -periodic solutions close to generating 
solutions with angular frequencies a;QjO)^ij> provided that H := luHq + eH\ + 
w 2 e$, the function Hq "projects" to the factor Mo, and \\Hq — -ffollc 2 + ll-^i — 
H\\\c 2 + ||^ — 'I'llc 2 ^ £ o- Moreover there exist at least 2 N ~ 2 such solutions 
counted with multiplicities. The phase orbits of these solutions are contained in 
some N -dimensional torus A that is O(oj)-close (and even 0(uj 2 )-close in the case 
(0,0) = 0) to the torus A° := U7=iU%n S 1 x {( J ij! > )} with respect to a 
C x -norm. The intersection points of these phase orbits with the transversal surface 
S (i.e. with the transversal section to the two-dimensional phase orbits of generating 
solutions on A ) coincide with critical points of the function ^l^ns w hcre ^> is the 
generating function of the perturbed succession map g T ^ ~ + ■ M — > M , 

1 := Jo + six (see definition 2.1 below). 

Similarly to remark 1.2 (A) one can show that, for any collection (a°,b°) of real 
numbers a°-,6°- under consideration and for small enough < uj <C 1, the period 
T can take an arbitrary value of the form T ^ 27rao(a°, h°)/uj, hence the quantity 
ui 2 T can be arbitrarily small. Thus any of the nondegeneracy conditions (30) and 
(31) can always be fulfilled. 

Theorem 2.2 (symmetric relatively-periodic solutions). Suppose that, 
under the hypothesis of theorem 2.1, each of three systems: the model system (24), 
the unperturbed system (27) and the perturbed system (33) is reversible. Then the 
perturbed system (33) admits exactly 2 N ~ 2 symmetric (T, a) -periodic solutions that 
are 0(ui)-close (and even O(oj 2 )-close in the case d{F°j)(0,0) = 0) to the generating 
symmetric solutions with the angular frequencies under consideration. Each of these 
2 N ~ 2 solutions smoothly depends on the pair of parameters (T,a). Moreover the 
function V^r is an even function in the collection cp of angular frequencies fij\j^, 
and the phase orbits of the symmetric (T, a) -periodic solutions pass through the 
points if of the torus A with the property if = —(p. 

Theorem 2.3 (stable periodic solution). Under the hypothesis of theorem 
2.1, suppose that all the numbers dI £ 3 (I°j, 0, 0) have the same sign, e.g. negative 
(respectively positive), and at least one of the following two conditions holds. The 
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first condition consists of the strong nondegeneracy condition: 



a $ [-duj 2 T, Ciw 2 T] + ttZ 



and of the following condition of having the same sign: all the signs 



rjij := sgn f^Tr 



d 2 g tJ (/°,o,o) 

9{qij,Pij) 2 



) 



1 $J i $J n, ^ j ' n. 



(34) 



are equal. The second condition is the following more delicate one: for any set of 
real numbers onj G R, 1 ^ i ^ n, ^ j ^ rij, such that 



the sum of any two, possibly coinciding, numbers of the set does not belong to the 
set 27rZ. Then the ((T, a) -periodic due to theorem 2.1) phase orbit of the perturbed 
system (33) passing through any nondegenerate critical point of local minimum 
(respectively maximum) of the function ^l^ns * s or bitally structurally stable in 
linear approximation (see definition 3.1). 

In all the next statements of the present section, the nondegeneracy condition 
(30) is not assumed to be fulfilled. 

Theorem 2.4 (necessity of the nondegeneracy condition). Suppose that, 
under the hypothesis of theorem 2.1, the number u) G (0,wo] an d ^ e collection of 
angular frequencies Qjj = fiyfi?-, 0, 0) satisfy all the conditions apart from the 
nondegeneracy condition. Suppose that a G 27rZ and that a smooth function Rq 
on M "projects" to the factor Mq. Let us fix the two-dimensional torus 7 C A° 
corresponding to a T -periodic solution of the model system (24). Suppose that some 
(and, hence, any) point of 7 is not a critical point of the function (Rq) that is 
obtained by averaging the function Rq := Ralf^^^ 1 (h (I° 00)) a ^ on 9 T-periodic 
solutions of the model system. Then, for any real number D > 0, there exist a real 
number fiQ > and a neighbourhood Uq of the projection p(j) of the two-dimensional 
torus 7 to Mq such that, for any e > 0, uje/fiQ $J fi ^ and \T — T\ + \a\ ^ Dfj,, 
the following holds. The neighbourhood U :— Uq x Mi 0/7 does not contain any 
(T, a) -periodic solution of the perturbed system (33), provided that Ho — Hq + /xi?o, 
the function Rq "projects" to the factor Mq, the function H\ "projects" to the factor 
Ah, \\R - R \\ C 2 no, and ||$|| C 2 < D. 

2.1. Scheme of the proof of theorems 2.1—2.4. Let us describe two main 
stages of the proof of theorems 2.1-2.4. 

Stage one is based on generalizing the averaging method (which was initially 
[12, 13, 7] introduced for Hamiltonian systems) to the case of 0-Hamiltonian systems 
with slow and fast variables. (A similar generalization see in [14].) It studies the 
unperturbed system (27) taking into account that it is close to a "super-integrable" 
model system (24). At first, one describes (T, a)-periodic solutions of the unperturbed 
system (27) that are close to the generating solutions of the model system (24). At 
second, one studies the linearization of the "succession map" for such solutions (see 
theorem 2.5 and corollary 2.1). 



am = Vio&, Wij - r)ij( a + A j:) cj 2 T)| < C 2 u 3 T, 1 < i ^ n, 1 < j < n. 
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As is well known, any symplectic linear operator decomposes into the direct 
product of "indecomposable" symplectic operators. Each of the latter operators is 
given by a "standard" symplectic matrix called a Jordan- Kronecker block in some 
canonical (i.e. symplectic) basis. 

Theorem 2.5 (Jordan-Kronecker blocks of a linearization of the 
UNPERTURBED SUCCESSION MAP) . Suppose that, under the hypothesis of theorem 
2.1, the number uj G (0,Wo] an d the set of angular frequencies = ( , 0, 0) 
satisfy all the conditions apart from the nondegeneracy conditions. Then there exists 
a N -dimensional torus A that is 0(iu)-close [and even O(oj 2 )-close in the case of 
d(-Fj°-)(0, 0) = 0) with respect to a C x -norm to the torus A° and is formed by the 
phase orbits of (T, a) -periodic solutions of the unperturbed system (27). Moreover, 
for any point x G A, there exist canonical frames eiji,eij2,eij3,eij4 in the tangent 
spaces T Xij Mij (where Xij :— pr^x), pr^ : M — > Mij is the projection) such that 
the linear part d(g^ Ho h 1 +lj 2 <s97o°'ii^ x ^ °f ^ e unperturbed "succession map" at the 
point X with respect to this frame is given by a blockwise lower-triangular matrix 
with the diagonal blocks 



( 1 J 1 




V o 



1 








cos a 






sin a 



'j 



\ 



J 



1 i ^ n, 
j ^ Ui 



Here are some real numbers such that = r\i^OL, \onj — r/ij(a + AijUj 2 T)\ ^ 
C2W 3 T for 1 ^ j ^ ni, while the numbers Ay and the signs r/ij G {1, —1} are the 
same as in (32) and (34). Furthermore all non-diagonal blocks vanish, apart from 
those blocks whose row and line correspond to the factors Mio and My (respectively) 
in the direct product M = n£=i IIj=o • ^he vectors eijk are bounded (k = 
1,2,3,4), and the relations e^i = d/difij, ey2 = d/dlij, e^ G R > od/dqij are 
either exact for j = or hold up to 0(uj) for 1 < j < Uj. 

In particular, onj ^ (mod 2tt) for any i,j if the nondegeneracy condition holds; 
aij ^ (mod 7r) if the strong nondegeneracy condition holds. 

Let us explain why the matrix in theorem 2.5 is a blockwise lower-triangular 
matrix. Unlike to the model system (24), which is a direct product, the unperturbed 
system (27) is only a semi-direct product. This is why, in theorem 2.5, the linear 
part of the unperturbed succession map g^ H ^ 1+a ,2$5/ Q / 1 is given by a blockwise 
lower-triangular matrix, but not by a blockwise diagonal matrix. Notice that the 
diagonal blocks in theorem 2.5 have the form exp(o; 3 TVij) where 



( 


^(/°-,o,o) 





o \ 

























V o 







o / 



1 ^ j ^ n, 

SC j ^ Hi 



(35) 



Q i0 = r)i tt l0 , - r)ij(Q.ij + A i3 uj 2 )\ < C 2 w 3 and a tj = O^T (mod 2tt) for 

l < j < m. 

From theorem 2.5, one can easily derive the following its refinement. 
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COROLLARY 2.1. Suppose that, under the hypothesis of theorem 2.5, the model 
and the unperturbed systems (24) and (27) are reversible. Then a ((T, a) -periodic 
by theorem 2.5) solution j(t) of the unperturbed system with 7(0) G A is symmetric 
if its initial point 7(0) =: {<Pij,Iij,qij,Pij}i,j satisfies the conditions ipij = —ifij 
(mod 271"). Hence the points x := 7(0) and x' :— g I ° l j 2 ('y(T/2)) are fixed under the 
involution J. Moreover, there exists a unique tangent frame e(x) = {e^i, ejj2, 6^-3, ejj4 
ai £/«e point x that satisfies the conditions of theorem 2.5 and consists of eigenvectors 
of the involution dJ{x) with eigenvalues {— 1, 1, 1, — l}ij. In the frames e{x) and 
e{x'), the operator A 1 = A'(x) = d{g^ o Hi+hj2 ^"(^){x) : T X M T X ,M is given 
by a blockwise lower-triangular matrix that is analogous to the matrix in theorem 
2.5 and has the diagonal blocks exp((cj 1 ~ : 'T/2)Vj 3 ), see (35). 

Stage two of the proof of theorems 2.1-2 .4 is based on generalizing the "method of 
generating function" (which was initially [4, 13, 15, 7, 16] introduced for Hamiltonian 
systems) to the case of an e-Hamiltonian ("perturbed") system that is C 2 -close to 
a O-Hamiltonian ("unperturbed") system. It studies T-periodic trajectories of the 
perturbed system in a neighbourhood of a "nondegenerate" compact submanifold 
A formed by the phase trajectories of T-periodic solutions of the unperturbed 
system. At first, one proves that the intersections points of T-periodic trajectories 
of the perturbed system with the "transversal surface" £ (see the formulation of 
theorem 2.1) coincide with critical points of the function ^| Sn x- At second, one 
studies the linearization of the "succession map" for these solutions. Here ^ is the 
generating function of the perturbed succession map (see definition 2.1), and A is a 
submanifold that is C 1 -close to the submanifold A. One proves the fact from above 
similarly to the case when both systems (the unperturbed and the perturbed ones) 
are Hamiltonian [13, 7]. In more detail, one constructs the submanifold A in the 
same way (since the construction does not use that the systems are Hamiltonian). 
Further, one proves that the submanifold A is nondegenerate and that critical points 
of the function ^^p^ coincide with critical points of the perturbed "succession map" 
via the results of the first stage. 

Definition 2.1 (generating function). Let e > and A : M — >■ M be 
a symplectic self-map of a symplectic manifold (M, ui) = (Mo x Mi,u>o + £<^i)- 
Denote Viji := ifij, x>ij% := qij ("coordinates"), Uiji :— lij, u%j2 '■— Pij ("impulses"). 
Denote by a the differential 1-form of the type of (u' — u)dv + (v — v')du' on M. 
More precisely, we define a by the formula 

n 2 / 

a ( x ) : =^2^2[ Kofc - M *ofc) dv mk + (v Mk - v' t0k ) du' m + 

i=l k=l \ 

rii 

+ £ Y1 ~ u ^ k ) dv ^ k + ( v ^ k ~ v 'vk) du 'ijk) 

where A : x = (v, u) 1— > A(x) =; (v', u'). In other words, 

a(x)£ := u(A(x) - x, dP(x)£), £ e T X M, (36) 

where P : x — (v, u) 1— » P{x) := (v, u'), and u)(^,n) denotes the value of the 
symplectic structure u> = YJi=i J2k=i( du i0k A dv iok + eJ2]Li du ijk A dv tjk ) on the 
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pair of vectors ^,tj e T X M . A function ^ = \&(x) will be called the generating 
function of the map A if it satisfies the condition 

d&(x) = a(x), x £ M. (37) 

Let us show that such a function ^> exists, i.e. the form a is exact. One easily 
shows that the integral of the form a along any closed curve equals the integral of 
the symplectic structure u) along some two-dimensional torus. The latter integral 
vanishes, since the symplectic structure under consideration is exact (being the 
standard symplectic structure on M = T*Q). Thus, the integral of the form a 
along any closed curve vanishes. This proves that the function ^ is well-defined up 
to an additive constant. 



§ 3. Normalizing the N + 1 body problem of 
the type of planetary system with satellites 

Consider the problem about the motion of a system of N + 1 particles attracting 
by Newton's law on a Euclidean plane E 2 , N ^ 2. The particles attract with the 
Newtonian potential U from (2). The equations of the motion have the form (1). 

The configuration manifold of the problem under consideration consists of all 
sets of radius vectors r^ <E E 2 , ^ i ^ N, such that r, ■ ^ rj, ^ i < j ^ N. In 
particular, the planar N + 1 body problem has 27V + 2 degrees of freedom. 

The system of equations of the motion in the cotangent bundle of the configuration 
manifold is a Hamiltonian system with the Hamiltonian function H being the total 
energy of the system: 

JV 2 

H = G + U where G = V #S U = - V (38) 



2lli 

i=0 ^ 0^i<j^N 



and the symplectic structure 



JV 

U = dp A dq = J2( d Pl A d ll + d Pf A d( ll)- 

i=0 

Here q = (r , . . . , r N ) = (rj, r 2 , . . . , r x N , r 2 N ) are coordinates and p = (p , . . . , pjy) = 
(Pq,Pq, . . ■ tPpf,p%) are impulses (i.e. the conjugate variables to the variables q), 
while G and U are the kinetic energy and the potential energy (2) of the system. 

3.1. Reducing to a system with 2N degrees of freedom. Consider the 
"reduced" problem, which is obtained by transferring to the coordinate system 
centred at the centre of masses 

c = c(r ,...,rjv) = — ^ 

L i=0 Mi 

of all particles. The configuration manifold Q of the obtained problem is an open 
subset of the 2 TV-dimensional vector space {(ro, rjv) £ (E 2 ) N+1 \ c(tq, . . . , rjv) = 
0}: 

Q = {{r ,...,r N )£ (E 2 ) N+1 \ c(r , . . . , r N ) = 0, r« ^ r^, i < j}. (39) 
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The variables q' = (ri, . . . , rjy) are taken as coordinates on Q. 

As a phase space of the problem, we take the cotangent bundle X :— T*Q to the 
configuration manifold. The phase space is 4A r -dimensional. Its natural coordinates 
are the variables q' (coordinates) and their conjugates p' (impulses). From now on, 
the accents will be omitted. 

3.2. Stating the problem. Our general problem is to find (T, a)-periodic 
solutions of the system described above, see definition 1.1. 

We observe that the motions with respect to the rotating coordinate system with 
angular velocity lj\ are described by the Hamiltonian system with the Hamiltonian 
function H — uj\I where 

N 

/:=X>i,Pi] (40) 

4=0 

is the "integral of areas", also called the angular momentum [17]. Here [q, p] is 
the oriented area of the parallelogram spanned by the vectors q and p. Thus the 
problem is equivalent to finding T-periodic solutions of the Hamiltonian system 
with the Hamiltonian function H — u\I where uji — a ~^ k , k being any integer. 
Due to (9) and (10), we can define uii to be the angular frequency of any of the 
planets or satellites, for example the angular frequency of the first planet. 

3.2.1. The parameters of the problem. Let us enumerate the particles of the 
system in a more convenient way (see (3) and (4)). Namely, we assume that one of 
the particles ro (the Sun) has mass /xo = 1, and the masses of the other particles (n 
planets and N — n satellites r^) have the form (3) where < /i, v -C 1 are small 
parameters and mi, mij > 0, 1 ^5 i ^ n, 1 ^ j ^ rii, are the parameters of the form 
(4). We will say that the ith planet and all its rii satellites form the ith satellite 
system. 

As above (see (8)), let us introduce "relative coordinates" in the configuration 
space Q (see (39)). Namely, we draw radius vectors 

Yij '■= r ij ~ r >> 1 < J < n U (41) 
from each planet to its satellites (if any). We draw also the "scaled" radius vector 
x, := (a - r )/R, 1 < i < n, (42) 
from the sun Tq to the centre of masses Cj of the ith satellite system where 



miY; 



3.2.2. Symmetric solutions. Let us show that conditions 1 and 2 in definition 
1.2 of a symmetric solution are equivalent. For this we will use the invariance of 
the total energy H of the system under the following two involutions Si and S in 
the phase space T*Q. 

Let I be any line in the plane of the motion going through the centre of masses 
of the system of particles. Define the following three transformations in the phase 
space T*Q preserving the total energy H of the system: 
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1) the canonical involution Si :T*Q — > T*Q corresponding to the axial symmetry, 
i.e. to the self-map of the configuration manifold Q sending all particles of the system 
to their images under the symmetry with respect to the line I; 

2) the anti-canonical involution S ("reversion of time") sending each pair (q, p) € 
T*Q to the pair (q, — p) where q and p are the sets of "coordinates" and "impulses" 
of all particles of the system; 

3) the anti-canonical involution J; := SSi = SiS. 

Each of these transformations is an involution, i.e. coincides with its inverse. 
The first involution is canonical, i.e. preserves the canonical symplectic structure 
dp A rfq on T*Q. The second and the third involutions are anti-canonical, i.e. 
they affect the symplectic structure by changing its sign to the opposite. Thus 
all three involutions move trajectories of the system to trajectories, moreover the 
first involution preserves the time on trajectories, while the second and the third 
involutions "reverse the time". 

A solution satisfies the first (respectively the second) condition of symmetry if 
and only if the point of the phase space corresponding to the time to (respectively 
any time t € R) of this solution is fixed (respectively is mapped to the point 
corresponding to the time 2i — t) under the involution J; = SSi = SiS. This shows 
the equivalence of the conditions 1 and 2 in definition of a symmetric solution. 

A solution j(t) is symmetric and (T, a)-periodic if and only if its points 7 (to) 
and j (to + T/2) are fixed under the involutions Ji and Jr o/2 (i) respectively where 
R v : R 2 — > R 2 denotes the rotation by the angle ip. 

3.2.3. A stable relatively-periodic solution. Suppose that a Hamiltonian system 
(X, u>, H) is S^-symmetric with respect to the Hamiltonian action of a circle S 1 on X 
via the Hamiltonian function I. Then the function / is a first integral of the system. 
Consider the flow g* H _ u j ■ X — > X, t e R, of the system with the Hamiltonian 
function H — The map A := g"ff_ u j will be called the succession map. and its 
linear part dA(x) at a fixed point x will be called the monodromy operator at this 
point. 

Let us define a "reduced" succession map for the two-dimensional torus 7 corresponding 
to a (T, wiT)-periodic solution. Let S C X be a small surface of codimension 2 that 
transversally intersects the torus 7 at some point xo = 70 E. Consider the restriction 
of the system to a regular common level set 

Xhj '■= {H = const, I = const} D 7 

of the first integrals H and /. Consider the small surface a := E n Xhj of 
codimension 2 in Xhj, which transversally intersects the torus 7 at the point 
xo = 7 n a. Consider the two-dimensional foliation on the manifold Xhj whose 
fibres are invariant under the (commuting) flows of the systems with Hamiltonian 
functions H and /; this condition uniquely determines fibres. Take the self-map A 
of the surface a sending any point x € a to the "next intersection point" of the fibre 
containing the point x with the surface a. In more detail, the map A : a' — > a is 
defined in a sufficiently small neighbourhood a' C a of the point x in a, it is "close" 
to the map A\ a i = g H -ij 1 i\<y' an d has the form A(x) — 9^-^ < x ) 1 Here /1 and 
fi are some smooth functions on a' defined by the conditions A(x) € er, x E a', 
fi(xo) = w i ; f2(xo) — T. The map A will be called the reduced succession map (or 
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the Poincare map), and its linear part A = dA(xo) at the point xq will be called 
the reduced monodromy operator corresponding to the torus 7. 

Recall that a linear operator A is called (Liapunov) stable if the norm of the 
operator A™ is bounded for n — > 00. As is well-known [17], the transversal surface 
(j is a symplectic submanifold and the self-map A of this surface is also symplectic. 
In particular, the reduced monodromy operator A is symplectic too. A symplectic 
operator A is called structurally stable if it is stable and any symplectic operator 
that is close enough to A is stable too. 

Definition 3.1. The two-dimensional torus 7 and the corresponding relatively- 
periodic solution will be called orbitally structurally stable in linear approximation 
(OSSL) (respectively orbitally stable in linear approximation on a common level 
surface of the first integrals of energy and angular momentum (OSLI)) if the 
reduced monodromy operator A = dA(xo) corresponding to the torus 7 is structurally 
stable (respectively stable). The torus 7 will be called isoenergetically nondegenerate 
(IN) if 1 does not belong to the spectrum of the reduced monodromy operator A, 
i.e. 1 ^ spec A. The torus 7 is called orbitally stable in linear approximation (OSL) 
if the linear operator B = dB(xo) is stable, where B : £' — > £ is the map defined 
similarly to the Poincare map A : a 1 — > a. 

Definition 3.2. An eigenvalue A e C of a symplectic operator A is called 
elliptic [17] if it satisfies one of the following equivalent conditions: 

1) the quadratic form Q(£) = a;(A£, £) is (positive or negative) definite on the 
maximal invariant subspace where the operator A has no eigenvalues apart from A 
and A; 

2) the Hermitian quadratic form ^) is (positive or negative) definite on 
the complex eigensubspace with eigenvalue A of the complexified space. 

The quadratic form Q is called the generating function of the symplectic operator 
A (see also definition 2.1). 

Proposition 3.1. (A) A symplectic operator is stable if and only if it is diagonalizable 
over C and all its eigenvalues belong to the unit circle in C. 

(B) A symplectic operator is structurally stable if and only if all its complex 
eigenvalues are elliptic. □ 

Let us mention some relations between the stability properties introduced above 
of an invariant two-dimensional torus 7: 

1) the following implications hold: IN <= OSSL => OSL =>■ OSLI. (The second 
implication is an important property of tori having the OSSL property; it follows 
from property 3 below. The first implication follows from proposition 3.1(B). The 
third implication is obvious. The inverse implications are in general false); 

2) if all eigenvalues of a symplectic operator A are pairwise different and lie on 
the unit circle in C then it is structurally stable, thus the torus 7 is OSSL; 

3) if the torus 7 is isoenergetically nondegenerate (IN) then it is included into 
a smooth two-parameter family of isoenergetically nondegenerate two-dimensional 
tori 7jj.j where parameters of the family are values of the first integrals H and /. 
If the torus 7 is OSSL (and, hence, IN) then all the invariant tori of this family are 
also OSSL. Hence OSSL =*> OSL (and not only OSLI). 
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We stress that, if the torus 7 is OSLI and IN, then the other tori of the family 
do not need to be OSLI, thus the torus 7 does not need to be OSL. 

Thus, among the stability properties introduced in definition 3.1 for a periodic 
solution, the strongest one is the OSSL property, which is studied in the present 
paper (see theorems 1.2(B) and 2.3). 

3.3. Describing the problem about planetary system with satellites 
via an e-Hamiltonian system. We define a smooth function F = F(x, y) on 
(R 2 \ {0}) x R 2 (or on (R 3 \ {0}) x R 3 ) by the formula 

F(x, y) = F 00 (x, y) := — - . (43) 

The function F will be called the Hill potential (or the "limit potential of the action 
of the Sun to a satellite"). As lemma 3.1 below shows, the unperturbed field of 
accelerations on the configuration manifold Q (which can be obtained in the limit 
when fi, v — > and uj > is fixed) is described by the following system of equations 
in the coordinates (41), (42): 

1 ^ i ^ n, 

(44) 

" w 2 ff (x^yy), l^j^m. 

For given i, j, the first equation of the system (44) describes the Kepler problem for 
the ith planet, while the second equation describes the so-called Hill problem for its 
jth satellite when the planet moves in a given way by virtue of the first equation. 

Remark 3.1. (A) The Hill potential F(x, y) is in fact the third coefficient of 
the power series of the function — | x _|T py | m the variable p at zero: 




5==:^-P%f-pXp(x,y), N.-o 



Py\ v/x 2 + 2p(x,y) + p 2 y 2 |x 

3x 2 y 2 -5(x,y 



^o,p(x, y) = F(x, y) - p{x, y) 



2|x|7 



_ p2 3xV-30x 2 (xy) 2 y 2 + 3 5 (x, y)4+ p-^Q. (46) 

8|x| 

A more general analytic potential 

Fgjx, y) - eF ,-e P (x, y) + (1 - 9)F , {1 _ e)p (x, y) (47) 

appears in the three-body problem (with < 6 < 1, fi, lj, p > and (8)) and in the 
restricted three-body problem (with 9 — 0, p,tu,p > and (8)), see §3.3.2 and (70). 
We remark that ff (x,y) = x2y ^^ y)x . 

(B) The unperturbed system (44) shows that the variables x.; and yy are automatically 
slow and fast variables respectively, provided that uj is small. 

Let x,, yij, rjij be the natural coordinates in the space X = T*Q corresponding 
to the coordinates x^, y^ in the configuration space Q, 1 ^ i ^ n, 1 ^ j ^ (see 
§3.1, 3.2.1). 
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Denote 



= nii + v } rriij, nii = — — , m tJ = — -, m l3 = f , (48) 

^— ' 1 + arm mi m l + i/ma 



the total mass of the ith satellite system, and the "reduced" masses of planets and 
satellites. Introduce the following functions on T*Q: the Hamiltonian functions 

ft €j "H = _ m,m, 3 

2rrii |xj| ' J 2mij \yij\ 

of the Kepler problems, the angular momenta 

IiO = [xi,€i], Iij = [yij,Vij], K K i, l^j^n* (50) 
of planets and satellites, and the "perturbation functions" 

1 ^ i < i' ^ n, 1 ^ j < j' ^ n,, of planetary system and satellite systems 
respectively (corresponding to pair- wise interactions of planets, respectively satellites 
of the same planet). 

As a "perturbation potential", let us consider the function 



n 



$ = $(x«,y««,m«,m«<.,/z, v,p) := y]mj®i + p mim v §a> (52) 

in the configuration variables Xj, € R 2 and parameters fhi, fhij,p, v, p £ R. Here 
the functions $i = $j(xj, y«, m», z/, p) and $jj/ = <&jj/(xj-Xj/, y«, y^,, mj*, m*/*, v, p) 
are defined by the formulae 

2^ 1 m l /m l ^\ mij/fhi 

l/p $i := -. r - i jrr _ v / 1 F7> ( 53 ) 

|xj| |xi-i/pdi| ^ |Xj + pyij - pvdi\ 

2 _ _ ffiifTii' ^ vm v mij 

vp rrumv^iv := , - > , ■ r ^— - 

|xj - Xj/| ^ |xj - Xj/ + pyy - ^(d, - 0j/)| 

mirrii' x ^ vmiini'ji 



E 



Here 



|x, - Xj/ - z/p(5j - 5j/)| ^ | - Xj + Xj/ + py»'j' + /w(5j - 5j/)| 

n, "«/ 2 

_ ^ m '3 m ''i' (54) 
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is the radius vector drawn from a planet to the centre of masses of the system of 
its satellites, multiplied by (E/=i mij)/ffii- One easily shows (see (45)) that the 
function $i is analytic in all its variables in the region 



1 + v E "Hi' + °> H 1 + M E 



ty' | 



3'=1 



i'=i 



while the function is analytic in all its variables in the region 



|p|(l + M£# =1 



|y«l < ¥ 



1 < i < rii, 

1 < / < ni. 



(56) 



(57) 



I |p|(l + klE^'=il^ 
The functions 3>^/ are expressed in terms of 3>i, .. . , $ n as follows: 

777- ' 

= -r 1 - $i (Xi - Xj/ + , y« , m« , i/, p) + 4^/ (x»/ - , y^* , m^* , v, p)+ 



v — — $j (xi - Xf/ - pyj/ j/ + i//5<5j/ ,y iif ,m it .,v,p). 

TTUi 

j'=l 



(58) 



We set $i := if rij = (i.e. the ith planet has no satellites), and we set := 
if rii = rii> = 0. If rii = 1 (i.e. the ith planet is a double planet, 9i := vmn/(mi + 
win)) then we have = e * ( . 1 ~ 0i ) Fg i;P (xi,ya) and 

= *i'(x,/ - Xi,yi/*,mi/»,i/,p) + — r -^-F Si „(xi - Xj/ + upS^, ya)+ 

vm v jm v 



%(l - 6i) V ^-F 9tiP {m - ya 

3 =1 



Wi' + vp8 v ,yn 



Lemma 3.1 (equivalence of the iV+1 body problem to an e-Hamiltonian 
System). Let Q be the 2N -dimensional vector space formed by all configurations 
of N + 1 particles with masses (3) and the centre of masses at the origin in a 
Euclidean plane. Define linear coordinates on Q to be the collection of radius vectors 
x i>yij '■ Q ~^ IK. 2 > 1 ^ i ^ n, 1 ^ j ^ rii ( see (41), (42)). There exists a collection 
of linear functions ^,17^ : Q* — > M 2 , 1 < i < n, 1 < j < (impulses) having 
the following properties. In the coordinates Xj, y^, rj^j on T*Q = Q x Q*, the 
canonical symplectic structure u) = dp A dq, the Hamiltonian function H of the 
N + 1 6ocft/ problem and the first integral I of angular momentum (see (38) and 
(40)) have the form 



H = -H, 



I 



P 



I 



(59) 



where 



uj = (jJq + eu>i, 



(60) 
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u> = d£ A dx = ^2 d£i A dx i: = dr] A dy = ^ ^ drj^ A dy r] , 

i—1 i—1 j—1 

n n I m \ 

H = ^K i+f , Yl & = E E^" E s ™ • ( 61 ) 

i=l l<i<i'<n i=l \i=l IC^i'^i I 

n n rii 

Jo = [x, £] = £ J.o, Jx = [y,»7] = EE J «- ( 62 ) 

i=i i=i i=i 

see (49), (50), (51). i/ere </ie small parameters < u;, e, /i, z^, p <C 1 are related by 
the conditions p — a; 2 / 3 /! 1 / 3 and e = u; 1 / 3 ^ 2 / 3 !/ = vp 2 /uj, p — -k. The "perturbation 
potential" $ = $(x», y„», to*, to**, A*i /?) ^as the form (52), is an analytic function 
on the direct product of the regions 

j 1 + v Y%Lx + 0, ] 

\ \P\ (! + IH E"'=i Ky'l) |y«| < mm(| Xi |, ±mhv #l |x* - x*/|) 



1 i n, and satisfies the condition 

( \ 



n rii 

$U=p=o = EE TO *. 
i=i j=i 



F(xi,yy) +/i^TOi/F(xi -Xj/,yy) 



(63) 



see (43) . In particular, H\ = $ = m i/ie case rij = of a planetary system without 
satellites. 

Remark 3.2. Lemma 3.1 implies equivalences of the following (e-)Hamiltoman 
systems for u, e, \i,v,p> 0: 

(T*Q,lj,H) £* (T*Q,u,H) S (T*Q,T*Q ,p;wo,u;i;JJ,JJi+a; 2 $) e 

where Qo is the configuration space of planets, and p : T*Q — » T*Qg is the 
projection. The third of these systems, called i/ie unperturbed system, is not only 
Hamiltonian, but also e-Hamiltonian (see (22)). Hence it naturally extends to any 
nonnegative values u>, e, fi, v, p of small parameters (despite of the fact that the 
symplectic structure degenerates if one of the parameters vanishes) . For the limiting 
values of the parameters uj > and fi = v = (and, hence, e = p = Q), the third 
system becomes a O-Hamiltonian system (T*Q, T*Qo,p; ujq, w±; loHo, Hi + u; 2 <f>)° 
called unperturbed where Hq := iJo|^=o> JJi : = H\\ v —o, $ := $| M=I / =p= o- The 
system (T*Q, T*Qo,p; ujq, wi; ujHq, Hi) , called the model system, is w 2 -close to 
the unperturbed system. It follows from lemma 3.1 that the unperturbed system 
indeed has the form (44) in the configuration space Q. 

Due to (45) and (58), the functions <&i, in (53), (54) have the following form 
for \v\ < v 0) \p\ ^ p : 

rii 

$i = v^Fo^pfa^i) + y^^r L F ,p(xi,y i j - vdi) = 
mi ^ mj 
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3 = 1 



vF(x l ,5 l )+0(vp), 



(64) 



rrii mi/ 
j=i i'=i 

-vF(y.i - -x.i>,8i) - ^F(x. ( / - x*,^/) + 0{vp). 
Hence the "perturbation potential" $ in (52) satisfies the condition 



5=EE" 



i=l j=l 



\ 



i'=l 



(65) 



-E 



mjF(xi, Ji) +mE miF ( Xi _ x i'^i 



V 



i =1 

Z 7^2 



/ 



Q{vp), v,p^0, 



which implies (63), see (46). 

Let us explain a geometric meaning of the Hill potential -F(x, y) when the «th 
planet has m > 1 satellites (1 ^ i ^ n). Consider the function $i = <&i(Xj,yj*, m**, z/, p) 
defined by the formula (53) and called the "potential of interaction of all satellites 
of the ith planet with the Sun". Due to (64), the function $i\ p =o equals a linear 
combination of the functions F(xi,yij), 1 ^ j ^ ni, and F(xj,5j). 

3.3.1. The Poincare transformation in the n + 1 body problem. In order to prove 
lemma 3.1, we will explicitly construct the variables of impulses £ it r}^ and will show 
that the function H in (59), (60), (61) equals the total energy of the system. One 
easily shows that those summands in H that do not depend on the impulses give 
the potential energy U. 

Let us compute the kinetic energy G. We will explore the fact that the transition 
from the coordinates (i"o, ri, . . . , r^r) in the configuration space to the coordinates 
Xj, yij (see §3.2.1) can be done by applying twice the following transformation called 
the Poincare transformation. 

Let us consider the configuration manifold Q of a planetary system (i.e. the 
system of n+1 particles in a Euclidean plane). It consists of all ordered collections of 
radius vectors r , ri, . . . , r„ with associated masses cq = 1, Ci = Ami, ■ ■ ■ ,c n = Xm n 
where < A 1 . Thus the manifold Q is naturally identified with the vector space 
R 2 (™+ 1 ) with coordinates ro, ri, . . . , r„. Consider the linear transformation L = 
Lc , Cl ,...,c n in R 2 ( n+1 - ) that corresponds to introducing the new linear coordinates 
on the space Q corresponding to the following collection of radius vectors: 



ro = 



r + ciri + 



C-nXr, 



l + Ci 



r x = ri - r , 



r n — r„ — ro 



where rn = c := r "t? ri ^'"T l ~ c " r " is the radius vector of the centre of masses of the 
system. 
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Definition 3.3. The transformation L = £ Co ,ci,...,e„ is called the Poincare trans- 
formation on the configuration manifold of the planetary system. 

Actually one could consider another transformation, namely the Jacobi transformation 
?o = c, f i = ri — c, . . . ,r„ = r„ — c. But this transformation would lead to more 
awkward formulae. Moreover it would bring us to a desired result only in the case 
of a usual planetary system, i.e. having no satellites. 

Consider the dual space Q* , i.e. the space of all linear functions on the space Q 
(or, equivalently, the cotangent space to Q at its any point). This space consists of 
all collections Po, Pi, • ■ • , Pn whose each item is a linear function on the plane, 
i.e. a co- vector. In fact, we can define the value of the linear function corresponding 
to such a collection on the configuration (ro, . . . , r„) g Q to be X)"=o(P i ' r *)- 
It is clear that the nondegenerate transformation L on R 2 ( n+1 ) != Q induces a 
linear transformation L* on (M 2 (" +1 ))* = Q* . Denote the image of the collection 
Po, pi, . . . , p n under the transformation L* by po, pi, . . . , p n . 

2 

Consider the real valued function G = Y^i=o 2c~ °^ kinetic energy on the space 
Q* . Besides we consider the function / = X^Lot 1 "*' P»] of angular momentum on the 
space T*Q. Finally consider the function of the total impulse P = Po + Pi + - • -+p« 
on Q* whose values belong to the space of co-vectors, i.e. of linear functions on the 
plane. 

Lemma 3.2. Under the Poincare transformation L = £ Co =i,ci,...,c„ on the configuration 
space Q of the n+l body problem, the functions G and P on Q* transform as follows: 

2 

(A) The kinetic energy G = X)"=o 2c~ ^ as ^ e f orm 

~2 n ~2 i 

G = 2 t+£lH< ii ' + - +& - )2 (66) 

where cq — 1 + c\ + . . . + c„ = 1 + \(ui\ + . . . + m n ). The expression (66) can be 
rewritten as follows: 

c = if +£§+ E »*> («o 

where c t := c Ci/(c + a t ) = c H /(l + c t ), 1 < i < n. 

(B) The total impulse P = po + pi + . . . + p Tl transforms to the impulse of the 
"heaviest" particle: 

P = Po- (68) 

The function of angular momentum I = X)™=o[ r *'P i ] on phase space X = 
T*Q is L-invariant: I = X^o^'P?]' 

PROOF. Items A and B directly follow by substituting into the functions G and 
P the following explicit formulae for the transformation L* of impulses: po = 

1+Cl+ 1 ...+ c „ Po ~ Pi Pn, P* = P* + 1+Cl +... +c „ Po, 1 < * < n. 

The invariance of the angular momentum follows from its invariance under the 
transformation on the space X — T*Q induced by any linear transformation = 
y~]. aijVj on the space Q with a nondegenerate matrix ||ay||. The latter holds, since 
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the transformation of impulses has the form p; = ^ fc bikPk where bikdij = 8jk, 
hence 



n n n 



YJ Yi ' Pi l = E E t r J ' Pk]aijbik = Y^i r j^PA = L u 

i=0 i=0 j=0 k=0 j=0 

Proof of lemma 3.1. Let us consider the case of a planetary system without 
satellites. Observe that the transition from the coordinates r to the coordinates 
x is a composition of the Poincare transformation L and the homothety = 
Rx.i, 1 ^ i ^ N. By setting pi = \f~^£,ii 1 ^5 i ^ N, po = 0, one obtains 
from (67) the desired expression for the kinetic energy G. In fact, GR equals 

zJi=i 2m~ + ^^2i^i<j^N\^i> ^j) ■ Hence, in the partial case of a planetary system 
without satellites, the function H in (59), (60), (61) indeed equals the total energy 
G + U of the system. The symplectic structure dp A dr = v / )j,Rd£ A e?x also has the 
desired form, since \J \iR = . 

In the general case of a planetary system with satellites, we observe that a 
transition from the radius vectors (ro, rj., . . . , rjy) to the coordinates x», ytj can be 
obtained via performing the following three transformations. At first, one should 
perform the Poincare transformation L to the whole system (see above). At second, 
one performs the transformation L to each satellite system r.^ , $J j ^ rij. Finally 
one performs the "scaling" homothety r^o = c» = Rx.i, 1 ^ i ^ n, Fy- = y^, 
1 ^ j ^ n-i- For the sake of simplicity, we will assume that all planets have 
satellites, i.e. all rii are positive. 

Step 1. The first performing of the Poincare transformation L to the initial 
configuration variables gives, due to (66), 

~2 N ~2 -. ~2 n -. 

G = £ + £ ft + 5® 1 + - + ~ mf " £ + S S< + jft + - + ^ 

i — 1 1—1 

Here = f£ + 2# = 2^7 + E"ii is the kinetic ener Sy of the ith 

satellite system, and = p^ + Y^jLi Pij ^ s its total impulse. 

Step 2. Put po = and perform separately the Poincare transformation L to 
each satellite system. As a result, we have from the formula (67) 

2/ifn< ^ 2//i/my fim, 1 ^,^ n . 

where mj, m.y are as in (48). By the formula (68), we have = p i . 
By using the previous step, we obtain 

^ I 9;;/m. 




-2 
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where m< = , , m '_ , 1 ^ i < n. 

S^ep 5. Now perform the following "scaling" of coordinates and impulses: 



Yi = Rx h fij = y i:j , Pi = J — ^ p y = fluffy (69) 



for 1 i $J n, 1 ^ j ^ n^. This gives: 




n / *2 «< „2 2 \ 

^Ehi^E^^ E + 1 E «..«. 



By taking into account that pvuR — e, we have the desired formula: 

n I pi n« 2 ^ \ ^ 

"**=E hi^Ei^ E <^> E 

i=l V ! j=l y 1 l^j<j'^"i / l^i<i'^n 



Due to (69), the symplectic structure has the form 



ijj = dp A dr = dp A dr = \J pR A dxj + e Yj (£77^ A dy^ 



since, recall, e — f^fl- This symplectic structure has the desired form (59), (60), 
(61), since \J pR = In a similar way, one proves the formulae for the first 

integral of angular momentum. 

This finishes the proof of the main lemma 3.1. □ 

3.3.2. Relations of the unperturbed system, the Hill problem, the three-body 
problem and the restricted three-body problem. Recall that the unperturbed system 
(see (44) and remark 3.2) is a O-Hamiltonian system consisting of n independent 
Kepler's problems (for each planet) and N — n Hill's problems (for each satellite of 
each planet). In particular, the Hill problem coincides with the unperturbed system 
for N = 2 and n — 1, i.e. it can be obtained from the three-body problem (with 
gravitational constant g — 1/p and the relation (8) between small parameters) in 
a limit when p, — > and v — > 0. As we will show below (see remark 3.3), the Hill 
problem can be obtained from the same 3-body problem in a limit when p, — > 0. 

Let us fix a real number to > 0, let the gravitational constant be equal to g = — , 
and study a limit of the Hamiltonian vector field in lemma 3.1 when p — > (or, 
equivalently, when the parameter p = -i = cj 2 / 3 /! 1 / 3 approaches zero). Here v is 
arbitrary (and unnecessarily v — > 0). 

Denote 8 := 8± — vmn/(mi + vmn), m :— rh\ = mi + vm\\. Since N — 2 and 
n = 1, the planet is a "double planet" (ni = 1). Hence $i(x, y) = et - 1 ~ e > Fg :P (x, y) 
and 

~ , 9(1-9)^ , 
<P(x,y) = mi$i(x,y) = m F<j iP (x,y). 

Here Fg tP (x,y) is the generalized Hill potential (or the "potential of the interaction 
of the satellite and the Sun") , which is determined for < 9 < 1 and p > by the 
relation 

1 (I 1-9 9 \ 

^<" (X ' Y) ■ 9(1 -0)? V R _ ]x~~9pT\ ~ \ X +(l-9)py\) ' (?0) 
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see (45), (46). It follows from lemma 3.1 and the relations e = vp 2 /u, fh\ = 
to, mi = 1+ " m , TOn = mi ^ mi = m0(l — 6)/i> that the 3-body problem under 
consideration with positive parameters 0<#<l,//>0, lj > and p > related 
by the condition (8) is equivalent to the Hamiltonian system with the Hamiltonian 
function 

~ / | 2 fhi\ ( rf toiTOh\ 2 ~ , 

H = uj\ t4 + e ^ r-r^ + w £$(x, y) = 

V2tox |x|; v 2toh |y| / V ,y) 

^ ton I f] mm\\ 



= u (i + fj,m)- — |—r ) + £\ m + w "^n^.p^y/ 

V 2m |x|/ V 2rri ii |y| 

and the symplectic structure u> — d£ A dx + ec??7 A c?y. Hence, it is described by the 
following system of equations: 



|f = + pm) + 0(1 - 0)p 2 ^(x,y)) 

^ = -TO^-c^(x,y). 



(71) 



The power series expansion (45) of the function r^q^l m the variable p at zero 
shows that the generalized Hill potential Fe jP (x, y) extends on the region |x| > |py|, 
^ 6 ^ 1 to an analytic function in all arguments and satisfies (47). Therefore 
Fgfl = Fqq = F. This implies that the 3-body problem (71) continuously and 
analytically extends on the region of nonnegative values of parameters ^ 6 ^ 1, 
/i^O, u > 0, p ^ related by the condition (8). 

One obtains the following problems from the 3-body problem (71) for limiting 
parameter values: the restricted three-body problem for 9 — 0, // > 0, u> > 0, p > 0, 
(8), and the Hill problem for ^ 6 ^ 1, uj > 0, [i = p = 0. More precisely, in a limit 
when u) > 0, p — » 0, (8) (and, hence, p — > 0), the three-body problem (71) uniformly 
in 9 G (0, 1) tends to the following system of equations in the configuration space: 

{ d 2 x _ „2 x . 
d'y y 2dF ( ^ ( 72 ) 

i.e. to the Hill problem coinciding with (44) for N = 2, n = nj = 1. 

Remark 3.3. (A) The Hill problem (72), in contrast to the three-body problem 
(with masses po = 1, pi, pn), does not depend on the parameter 9 = an d, 
hence, on the ratio of masses of the planet and the satellite. The Hill problem 
(72) formally depends on the parameters to > and u> > 0. But, after the change 

~ 2/3 

t = ojt, y = ^1/3 y> x = x, it transforms into the similar problem with to = to = 1. 
The latter problem does not depend on any parameters and it is usually regarded 
as the Hill problem. 

(B) The Hill problem was initially [1, 3] obtained as a limit when p — > and 
(8) from the restricted three-body problem (see (71) for 9 = 0), but not from the 
three-body problem (71) itself. As we showed above, the three-body problem (71) 
tends to the Hill problem (72) when p — > and (8), provided that 9 £ [0, 1] is 
arbitrary (but not only 9 = 0). This was observed already by Brown [18]. 
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3.4. Normalizing the Kepler problem near circular orbits. It follows 
from lemma 3.1 and remark 3.2 that, forO<(i,i/«iij<l 1 the N + l body problem 
of the type of planetary system with satellites is w 2 -close the the corresponding 
model system (MqXM\, U)o+u>i,u)Hq+Hi), which is the direct product of independent 
Kepler's problems 

(Mio,^ A dXi,u}H i0 := uK^ , (.\ /,,.<///,, A dy lJ: H l3 := sf^j , (73) 

1 < i < n, 1 < j < nj, with the Hamiltonian functions (49). Let us show that 
the Hamiltonian function of any planar Kepler's problem has the form H+j — 
Hij(Iij,qij,pij) and satisfies the conditions (28) of theorem 2.1, with respect to 
some canonical coordinates fij,Iij,qij,Pij on (Mi 3 ,u,j). 

The planar Kepler problem is given by the Hamiltonian system 

M = T*(R 2 \ {0}), u = dp A dq, H = - ^ = G + u] . (74) 

2m |q| J 

2 , 

Here G — ^ and U = — t^t are the kinetic and potential energies of the system, 
q € K 2 \ {0} and p £ M 2 are the radius vector and the phase impulse of the 
particle, m, k > are parameters. The solutions of the Kepler problem having 
negative energy levels H are periodic. Since the Kepler problem is invariant under 
all rotations of the plane, it has the first integral of angular momentum / = [q, p] . 
Note simple properties of circular motions in the Kepler problem: 

1) For any r > 0, there is a unique (up to changing the direction of rotation) 
circular motion of the particle satisfying the system (74) along a circle of radius r. 

The angular velocity of this motion equals ft = ±y^-, while the energy H and the 

angular momentum / equal H = — 4p and I = mflr 2 = ±mVkr respectively. In 
particular, the values H and / depend monotonically on the value Q (for Q > 
or £1 < 0) and take all values in the domains H < and 1^0 respectively. We 
will assume that the parameters r, Q, H, I of a circular motion are related by the 
formulae from above. 

2) Circular motions correspond to the equilibrium (i.e. stationary) positions of 
the particle with respect to a rotating coordinate system with angular velocity f2. 
Therefore, for any a ^ 0, the solution of the Kepler problem corresponding to a 
circular motion is (T, a)-periodic with T = In other words, such a solution is 
^-periodic with respect to a rotating coordinate system with angular velocity O (as 
well as any angular velocity of the form (1 + ^■)fl where I is an integer). 

The Kepler problem of our interest, i .\/, ,. <///, , A dyy,ily-) for 1 < j ^ m, has 
the form (74) with M = M ij7 H = H ij: I = Jy, 

(q,P) = (YijiVij)) m^rhij, k = rrii + urriij, r = r^-, f 2 = fly. (75) 

(Similarly, the Kepler problem (Mjo, d^AdXi, H^) has the form (74) with M = Mjo, 
H = H i0 , I = ho, (q, p) = (xi,£j), m = m,, k = 1 + /zmj, r = Rt/R = pRi = 
(1 + jiimj^lQjol -2 / 3 and Q, = fl i0 = Ui/uj.) 

Lemma 3.3 (normalizing the Kepler problem). Let H and I be the Hamiltonian 
function and the first integral of angular momentum of the planar Kepler problem 
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(74). In the region {r > 0, p^ ^ 0} of the phase space of the problem, consider the 
coordinates ip mod 2ir, I, q, p of the form 

2rp r km 2 r 
ip = 1p , I = ffy, q = In — — , p = rp r - (76) 

pv> pip 

Here r, ip are the polar coordinates in the plane of motion, p r , p^ are the corresponding 
impulses. Then: 

(a) the coordinates (76) are canonical, i.e. U) = dl A dip + dp A dq; 

(b) in the coordinates (76), the Hamiltonian function H of the Kepler problem 
does not depend on the angular coordinate ip mod 2tt and has the form 



H I 2 +p 2 1 1 q 2 p 2 



+ o(q 2 +p 2 ) 



k 2 m 3 2I 4 e 2 i Pel 2I 2 2I 2 2I 4 
as q 2 + p 2 — > 0. Furthermore the involutions Si, S in §3.2.2 have the form 

Si(<p,I,q,p) = (-<p,-I,q,p), S(<p,I,q,p) = (<p,-I,q,-p). (77) 

One proves lemma 3.3 in a direct way. □ 
The canonical coordinates ip mod 2tt, I, q,p in lemma 3.3, as follows from their 
construction, are quite "similar" to the canonical coordinates -0mod27r,p^ = /, r,p r 
corresponding to the polar coordinates ip mod 2ir, r in the configuration space of the 
planar Kepler problem. For example, the involutions Si, S have the form (77) in 
the coordinates ip,p^,r,p r too: 

Si(ip,Pi,,r,p r ) = (-ip,-p^,r,p r ), S(ip,p^,r,p r ) = (ip, -pj,,r, -p r ). 

For any ft ^ 0, denote by 7n the phase trajectory of the Kepler problem 
corresponding to the circular motion with angular velocity Q. The invariant two-dimensional 
surface Un^o7fi m the phase space formed by all these trajectories will be called the 
surface of circular motions. This surface is smooth and consists of two connected 
components, each of which is diffeomorphic to a punctured plane and bijectively 
projects onto the configuration manifold under the canonical projection. 

We recall that the system describing the motion with respect to a rotating 
coordinate system with angular velocity £1 is a Hamiltonian system with the Hamiltonian 
function H — 0,1 . Thus the circular motions correspond to the stationary points of 
such systems. 

COROLLARY 3.1. The surface of circular motions of the planar Kepler problem 
is a region in the symplectic coordinate cylinder: 

|J 7 fi = {(</? mod 2^, 7,0,0) I 

neR\{o} 

For any Q G R \ {0}, the quadratic part (i.e. the quadratic form whose matrix is 
formed by the second partial derivatives) of the function H — £11 at any point of the 
circle 7n has diagonal form with respect to the coordinates (76): 

S 2 (H - ni)\ yn = ^5I 2 + n (iSq 2 + %-) . (78) 

mr z V II 
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§4. Deriving theorems 1.1—1.3 from theorems 2.1-2.4 

Due to lemma 3.1 and remark 3.2, the N + 1 body problem of the type of 
planetary system with satellites is equivalent to the e-Hamiltonian system (33) 
with small parameters < uj, e, p, v, p <C 1 related by the conditions e = w 1 / 3 'p 2 / 3 v 
and p = w 2 / 3 ^ 1 / 3 . Moreover the functions Hq, Hi, $ are S 1 -invariant, the function 
H = H + pR "projects" to M := T*Q , the function H 1 = H x + vR x "projects" 
to Mi :— T*Q\, and their "principal parts" equal the sums Hq = 53i=i-^»o an d 
Hi = Ya=i 2^=1 Hij- Furthermore, due to lemma 3.3, each summand has the 
form Hij = Hij(Iij,qij,pij) and satisfies the conditions (28). The perturbation 
potential $ is an analytic function in a neighbourhood of any torus A°, provided 
that the collection of angular frequencies fly satisfies the conditions (6) and (7) of 
"lack of collisions". The principal part $ := ^\^— u — p — of the perturbation potential 
has the form (26). 

So, the N + 1 body problem of the type of planetary system with satellites 
considered in theorem 1.1 is equivalent to an e-Hamiltonian system belonging to 
the class of "perturbed" systems in theorem 2.1. 

PROOF of theorems 1.1 and 1.2. Step 1. In theorems 1.1 and 2.1, the 
"relative resonance" conditions (9) and (29) on the collection of frequencies are 
equivalent. The nondegeneracy condition (11) from theorem 1.1 is equivalent to the 
nondegeneracy condition (30) from theorem 2.1. 

Let us suppose that the more delicate nondegeneracy condition (12) from theorem 
1.1 holds. Let us prove the nondegeneracy conditions (31) from theorem 2.1. 
The first condition in (31) is equivalent to the first condition in (12). In order 
to prove the second condition in (31), let us evaluate the number Ay in (32). 
Due to (63), we have i*y — .Fy(xj,yy) = mijF(Xi,yij). By construction, = 
(Fij(Xi, ■))\ H - 1 [H ij (i9.,o,o)) wnere x i = const, |xj| = Ri/R = \fl i0 \^ 2/3 (see theorem 
2.1). Let (F°j) = (F°j)(qij,pij) be the function obtained by averaging the Hill 
potential F°j — F°j (qij , pij) along the ^--periodic solutions of the Kepler problem 
(Mij , u)ij , Hij) for the satellite. By an easy calculation, taking into account (43) 
and corollary 3.1, we find the differential and the Hesse matrix of the function (F°j) 
at the point (0, 0): 

d(F o \ fn Q) n ^°)(0,0) _ o2 ( - 2 9/8 
^ )(0 ' 0) -°' d( qij , Pij ? ~ n ^[ 25/(84) 

This and (78) imply that 



Ar = ^Tr [ [ d 2g ^°' M d 2 (f;°)(0,0 ) 
n 2 _ Tr (( I \f -29/8 ^ _ <>7„ 



2fiy" VV % J V 25/(8/ 2 ) J J 4^ 
Therefore the second desired condition in (31) has the form 



a + A^T = a- ^-T # [-C 2 uj 3 T, C 2 lj 3 T] + 2ttZ, 

4:1 L-lj 
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i.e. it is equivalent to the second condition in (12) with the constant C := C*2- 

Thus, all the conditions of theorem 2.1 are fulfilled. Hence this theorem implies 
theorem 1.2(A). 

Step 2. Let us check that, for the N + 1 body problem under consideration, 
the model problem, the unperturbed and the perturbed problems are reversible. 
The construction of the functions H, Hq, Hi shows that they (and, hence, also <&) 
are invariant under each of the involutions Si and S from §3.2.2. Hence they are 
invariant under the composition J — SiS — SSi. Due to (77), the involution J = 
SiS = SSi acts component-wise in the form (y>y, Zy, <Zy,f>y) H> (— </3y, hj, Qij, ~Pij)- 
From here, by taking into account the J-invariance of the functions Hq, Hi, we 
obtain the reversibility of the model system, of the unperturbed and the perturbed 
systems. Hence theorem 2.2 implies theorem 1.1. 

Step 3. Let us derive theorem 1.2(B) from theorem 2.3. By lemma 3.3 or (78), 

d 2 H 

all numbers dji ' J (i?-, 0, 0) are negative. The sign rjij from (34) equals 

/ 3 2 tfy(I?,,0,0)\ 
Vij = sgn fiyTr — -- = sgnfty. 

Hence, by the first property of having fixed sign in theorem 1.2, all the signs r^o = 
sgnilio are the same, moreover rjij Ay < for 1 ^ j ' ^ rii. Suppose that the 
conditions (15) and (16) hold for C := C 2 . Then a $ 7rZ and 



2 ' 2 



vrZ, (79) 



T <£ [-C 2 uj 3 T, C 2 u 3 T] + ttZ (80) 

for 1 ^ j ^ rii and 1 ^ j' ^ . Consider any collection of real numbers ay , 
1 < i < n, < i < rij, such that 

= Vio a i \ a ij - Vij( a + Av,w 2 T)| ^ C 2 uj 3 T, H K n, 1 s$ j < rii. 

Then: 

1) the sum a^o + a^o = 2?yioa does not belong to 27rZ, since a ^ 7rZ; 

2) for 1 s$ j ^ m, the sum a^o+ay e {r)w+Vij) a +Vij ■ Ai J uj 2 T+[-C 2 u 3 T 1 C 2 w 3 T] 
does not belong to 27rZ because of (79); 

3) for 1 ^ j ^ rii and 1 ^ j' ^ n,-/, the sum ay + ai>j' G (r?y +r]i>j/)a + {rjij Ay + 
ri Vr & vr )uj 2 T + [-2C 2 uj 3 T, 2C 2 ui 3 T] does not belong to 2ttZ because of (80). 

Thus, the sum of any two, possibly coinciding, numbers of the set ay does not 
belong to the set 27rZ. Hence, the hypothesis of theorem 2.3 holds, and therefore 
this theorem implies theorem 1.2(B). □ 

4.1. Necessity of the nondegeneracy condition a 7^ 0. In this section, we 
give a more exact definition of the notion "almost any" in theorem 1.3 (meaning 
"any" for N ^ n = 2, see corollary 1.1$)) and of the subsets M C M sym C R%. 
Besides we derive theorem 1.3 and corollary 1.1$) from theorem 2.4. 
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Consider the sequence of positive real numbers 

c K := — - / , K /i 2 \ 3/2 cos(nt)dt where r K := ( ) , (81) 

7T J-tt (r« + ^ - 2cos<) J / 2 \ / 

k 6 2 \ {-1,0}. It will interest us only up to a nonzero multiplicative factor. Here 
c K > due to the following properties of the integrand g K (t) in (81): it is an even 
function, g K (t) > 2\g K (2t + ^)\ for < t < g^, and fljt (t) > | 5ft (i + ^)| for 

^Tpf < ^ < ^ ^Utl ' The following property of the sequence c K is 

convenient to use for approximate computations: 

c K = C\k 2 + -(Ci + C 2 )k + o(k) ~ Cik 2 , |«| -> oo, 

where 

2 f 00 costrfi 1 f 00 costcft 

Cl " - loo (§T^ > ' C2 ^ * loo (fT^ > 

On the basis of numeric computations, which were carried out by A. B. Kudryavtsev 
in the interval \k\ ^ 1000, one can make a conjecture that the sequence c k /k 2 
decreases for < k — > oo and increases for — 1 > n — > — oo to its limit C\ s» 2.18. 
(This is true in the indicated interval.) Approximate values of c k /k 2 for small \k\ 
are given in the following table: 



K 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


C K /K 2 


17.55 


7.21 


5.04 


4.15 


3.67 


3.37 


3.17 


3.03 


2.92 


2.84 


K 


-2 


-3 


-4 


-5 


-6 


-7 


-8 


-9 


-10 


-11 


c k /k 2 


0.37 


0.72 


0.97 


1.16 


1.29 


1.39 


1.48 


1.54 


1.60 


1.64 



Put 

Kw := — , J^i', 1 < < n. (82) 

Then ^ = , Km = Kw + 1, thus the numbers K^i and Kw are either both 

integer or both non-integer, moreover k := Ku> G R\ {0,-1}. The number 
equals the ratio ^ of the angular frequencies of two planets along the circular 

orbits. Hence, due to Kepler's second law, the number r K = ( £ ^-) 2 ^ 3 in (81) equals 
the ratio of the radii of these orbits. Let us define the number c K G M for any eel 
as follows: either by the formula (81) if K € Z \ { — 1, 0}, or by the formula 

c K :=0 if k G K \ (Z\ {—1,0}). (83) 

Consider the following collection of complex-valued functions on the n-dimensional 
torus (S 1 ) n with angular coordinates (p = (cpi, . . . , ip n ): 

fw(cp)~K U >c Kll ,e iK «'^'-v>\ l^l', m,I'<n, (84) 

where i = >/— 1 G C is the imaginary unit and c K is defined in (81) and (83). 
The functions (84) are 27r-periodic in each argument, moreover they are equivariant: 

fll'ifl +t, <p n +t) = fw(<pi, ■ • .,<p n ), t£t, 

fll'((pi +U!it,...,tp n +U) n t) = e Wlt fu>((p-L,...,(p n ), t G R. 
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Definition 4.1. Let us fix the angular frequencies u>i of the planets satisfying 
the properties (9), (10), (6), (7). A collection of planets' masses [i{mi, . . . , m n ) £ 
R™q will be called unclosing for a phase point cp = {cp\, . . . , cp n ) € {S 1 )"", or simply 
cp -undo sing 1 , if at least one of the following numbers does not vanish: 

n 

fl(<p;mi,...,m n ) :=^2mi>fw((p) eC, U K n, (85) 
i'=i 

see (81)-(84). A phase point cp € (S 1 ) n will be called symmetric if it is fixed 
under the involution (5 1 )" — > {S 1 ) 11 , cp H> — cp, of the n-dimensional torus, i.e. 
its coordinates have the form cpi £ {0, it} mod 2ir, 1 ^ I ^ n. Denote by A4 
(respectively Ai sym ) the set of all collections of planets' masses fj,{mi, . . . , m n ) £ 
K™ that are i^-unclosing for any (respectively for any symmetric) phase point 
cp £ {S v ) n . 

Remark 4.1. The number of symmetric phase points equals 2™. For any phase 
point ip £ (S 1 ) n , the set of yj-closing collections of planets' masses is the intersection 
of a linear subspace of R" with M™ - Hence M sym is open in R™ . Moreover it 
is dense in R™ whenever it is nonempty. The subset M. C R> is open in R™ , 
since the torus (5 1 )™ is compact and the functions f\ = fi{cp; mi, . . . , m n ) are 
continuous, see (85). Suppose that Kw £ Z for some i ^ i', 1 ^ i, i' ^ n. Then the 
set of tp-closing collections of planets' masses is contained in a plane of codimension 
^ 2 (since the system of functions (85) is linear in (mi, . . . , m n ) and its rank is at 
least 2). Moreover any collection of masses with \ku> \c K ..,mii > ^ \Ka\c Kil mi 

is (^-unclosing for any phase point cp £ (S 1 ) 71 (i.e. belongs to A4). Hence the open 
sets M. C A4 sym are nonempty, thus "almost any" collection of planets' masses (see 
the paragraph before theorem 1.3) belongs to A4 sym . 

Consider the natural angular coordinates on the torus A°: 

x i y {cpu <fij, 1 < I < n, 1 < j ^ n t }, x £ A°. 
The following statement generalizes theorem 1.3. 

Proposition 4.1. Consider the N + 1 body problem of the type of planetary 
system with {or without) satellites, N ^ n ^ 2. Under the hypothesis of theorem 
1.1, fix the angular frequencies 0J\, ...,u„ of planets having the form (9), (10), (6), 
(7). Suppose that there exists at least one pair of planets with indices i ^ i' , whose 
frequencies are in a special resonance (21). In this case, one automatically has 
a = 0, Ha* £ Z \ {0,-1} and c K .., > 0. Fix the two-dimensional torus 7 C A° 
corresponding to a T -periodic solution of the model system. Let us suppose that the 
collection of planets' masses ft(mi, . . . , m n ) £ R" is cp -unclosing for some {and, 
hence, any) point x = {cpi, cpij} £ 7, see (85). 

Then, for any real number D > 0, there exist numbers /io, vq > and a neighbourhood 
Uq of the projection of the two-dimensional torus 7 to the phase space of planets 

X A collection of masses /i(mi, . . . ,m n ) £ R>o ' s cp-closing if and only if, for any index I = 
1, ... ,7ig the planar polygonal line An(ip) . . . Ai n (cp) C C is closed, provided that the segments of 

this polygonal line have the form A[[i (cp)A l ii+\{cp) = my hi' (f) f° r 1 $ J' ^ J — 1 and the form 

M,l'-x{<P)A w {cp)=mitf w {tp) for 1 + 1 sC I' ^ n. 
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such that the following property holds. For any values \x,v such that < (^) 3 ^ 
fi ^ fio (respectively < /! ^5 /j,q if there are no satellites), the direct product 
U of the neighbourhood Uq and the phase space of satellites does not contain any 
(T, a) -periodic trajectory of the N + 1 body problem under consideration, provided 
that the parameters T , a have the form 

|T-T| + |5| s^D/i. 

In particular, the assertions of theorem 1.3 hold. 

Let us show that theorem 2.4 implies corollary l.l(^) about a planetary system 
with two planets without satellites (N = n = 2). By lemma 3.1, the Hamiltonian 
function and the symplectic structure of the perturbed problem are 



ujH a = uj \ Ki + K 2 + nKi2) , <^o = d£i A dxi + d£ 2 A dx 2 . 

Here Ki — t^t is the Hamiltonian function of the Kepler problem corresponding 

to the ith pknet,^ = j^-, i = 1,2, K 12 = - 

Suppose that the ratio of frequencies ui\, ui 2 is rational, i.e. they have the form 

27rfci 2izk 2 



T T 

where ki, k 2 are nonzero integers, k\ ^ ±k 2 , T > 0. Then the solutions of the 
unperturbed problem (/i = 0) corresponding to the independent circular motions of 
planets with angular frequencies uj\, u 2 are T-periodic. Conversely, if the solution 
is T-periodic then the pair of angular frequencies (5) is proportional to a pair of 
integers, with coefficient Sr. 

Proof of corollary l.l(^). Step 1. Recall that the unperturbed problem 
(corresponding to fj, = 0) splits into two independent planar Kepler's problem. 
Hence its T-periodic phase trajectories form the six-dimensional submanifold 

Q = {Ki = const, K 2 = const} 

in the eight-dimensional phase space. In fact, due to periodicity of solutions of the 
Kepler problem with negative energy levels, the period of any its closed trajectory 
is a smooth (and strictly monotone, see §3.4) function in the value of energy. 

Let us find the averaged perturbation (Rq) , i.e. the function obtained by averaging 
the perturbation function Rq = Rq\q = (Ki 2 + ^£ )|e along the periodic solutions of 
the unperturbed problem. In more detail, let us find the differential of the function 
(Rq) at any point of the torus A° C 8. 

Step 2. With respect to polar coordinates ip, r on the plane of motion, we have 

~ Pr t +PlJ r ! fhi . 1 „ 

&i = 7P= > 1 = 

2rm n 

r\r 2 
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r 2 ri J yVf +r\- 2nr 2 cos(V>i - ^2) 

Step 3. Let us transfer to the coordinates tfi,Ii,qi,Pi, i — 1,2, from lemma 3.3. 
By this lemma, 

A° = { pi = p 2 = qi = 92 = 0, h = const, I 2 = const }, 

moreover the restriction of the linearized unperturbed system to the tangent bundle 
Ta<>0 = {£ I <£ii(£) = dl 2 (0 = 0} to 6 has the following form (with respect to the 
coordinates ipi, dtpi, dqt, dpi, i — 1,2, on this bundle): 

difi d(d<pi) d(dqi) dpi 1 d(dp. t ) 

— = {) > — = j— = -^ (86) 

i = 1, 2. Hence, at each point (<^i, y> 2 ) of the torus A°, the differential of the function 
i4fi2|e has the following form: 

Airs 1 \< \ ™-\™-2r\ri ( . ( Jpi dp 2 
d(Ki2\e){<Pi,<P2) = 5 sm^ia dip 12 + 2— 2— — 

n 2 v V h h 

— - cos tp 12 ) dqi + ( — — cos ^i2 ) dq 2 
r 2 J \n J 

^y- (cos tp r2 (dqi +dq 2 ) + sin</?i2 (d<pi 2 + -p- - ) j (S?) 



where ^12 := <p\ — (f 2 , T\ 2 '■= V r i + r 2 — 2rir 2 cos <^i 2 . The perturbation function 
has the form i?Q = Rq\& — (if 12 + |£ )|e- One easily shows that the contribution 
of the summand 7r£ 2 |e = \Y^i=i(Pri + Pv> i / r f))l© to the averaged perturbation 
(i?o) has a trivial differential at any point (ipi,ip 2 ) G A°, i.e. d(Rg)((pi,ip 2 ) — 
d(Kx 2 \@)(<pi,ip 2 ). 

Step 4- Consider the rational number n — n\ 2 = fc *^ fc — , see (82). If 

k ; G Z then define the number c K by the formula (81). 

By integrating the values of the co- vector d(Ki 2 \@) (see (87)) on the solutions of 
the linearized system (86), one obtains the following: 

1) if K G Q \ Z then the differential of the function (Rq) vanishes at any point of 
the two-dimensional torus A°; 

2) if k G Z then this differential has the following form at any point (ipi, ip 2 ) G A : 

d(R%)(ipi,ip 2 ) = mythi ( — ( cos(K^i 2 )d<7i - sdn^u) ^ 



(« + 1)C-k-1 / // . ,s s, - / / 1 -1 \ \ d P2 

cos((k + l)tpi 2 )dq 2 - sm k + l)^i 2 ) — 

r 2 \ h 

and, hence, it does not vanish (since c K ^ 0, c_ K _i 7^ 0, see (81)). 

Since any point of the torus A° is noncritical for the function (Rq), theorem 2.4 
implies corollary l.l(^). □ 
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PROOF of proposition 4.1. By lemma 3.1, the summands H , w in the decom- 
positions (60) of the Hamiltonian function and the symplectic structure of the N+l 
body problem have the form 

n n 

Ho = ^ Ki + n ^2 K u' > w o = ^2 d €i A dXl - 

i—1 l^i<z'^n i—1 

Here Ki is the Hamiltonian function of the Kepler problem of the ith planet that 
is similar to the Hamiltonian function K\, and is the function similar to K12. 
Put K «' = Zjf^T,' 1 < *,»' < n, see (82). 

From the equality (88) in the case N = n = 2, we immediately obtain the 
analogous formula in general case N ^ n 2: 

d{Ro)\ A ° =d(Ro)(<p 1 ,...,ip n )= Zii'(v>u<Pi')- 

l^i<i' 

Here £,* i i(<fi, ( fii') is the co- vector analogous to the co- vector (88), which is denoted 
by £12(^1 j ^2)- In more detail, we have 

n _ n ✓ 
d(i?o)U° = ~ ™"i' K ii' c Kui ( COs(Kiitiput)dqi 

1=1 i'=l 

where y^/ := cpi — ipii. This co- vector vanishes at those points of the torus A where 
the functions //, 1 ^ I ^ n, simultaneously vanish, see (85). Hence theorem 2.4 
implies the absence of (T, S)-periodic solutions in some neighbourhood of the torus 
7, provided that the parameters LJ,fx,e> are small enough and are related by the 
inequalities uje/fiQ /i /io- Due to the relation e these inequalities 

have the form w 4 / 3 // 2 / 3 ^/^ ^ M ^ i- e - the form oj 4 (v / /io) 3 ^ /1 ^ /xq. Therefore 
theorem 2.4 indeed implies proposition 4.1. □ 

Theorem 1.3 obviously follows from proposition 4.1, definition 4.1 of the subsets 
M C M sym C R™ , and remark 4.1. 
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